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Abstract 



We compute M = 4 Super Yang Mills planar amplitudes at strong coupling by consid- 
ering minimal surfaces in AdS 5 space. The surfaces end on a null polygonal contour at the 
boundary of AdS. We show how to compute the area of the surfaces as a function of the 
conformal cross ratios characterizing the polygon at the boundary. We reduce the problem 
to a simple set of functional equations for the cross ratios as functions of the spectral pa- 
rameter. These equations have the form of Thermodynamic Bethe Ansatz equations. The 
area is the free energy of the TBA system. We consider any number of gluons and in any 
kinematic configuration. 
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1 Introduction 

In this paper we consider minimal area surfaces in AdS space that end on a null polygonal 
contour at the boundary of AdS. Our goal is to compute the area of the surfaces as a 
function of the shape of the contour. Our solution to the problem consists of a system of 
integral equations of the thermodynamic Bethe ansatz (TBA) form [1] . The area is given by 
the TBA free energy of the system. 

Our motivation for this investigation is the study of scattering amplitudes in M = 4 super 
Yang Mills. Planar M = 4 Super Yang Mills is an integrable theory [2]. This means that 
if one finds the appropriate trick, one is going to be able to perform computations for all 
values of the 't Hooft coupling A [3, 4]. Finding the appropriate trick is usually tricky. Via 
the AdS /C FT correspondence, this problem amounts to solving the quantum sigma model 
describing strings in AdS*, x S 5 . The classical limit of this theory is simpler to analyze. This 
is what we do in this paper. We consider classical solutions for strings moving in AdS$. 
In the classical limit we can forget about the worldsheet fermions and the five sphere and 
study strings that are in AdS§. We think that the knowledge of these classical solutions will 
be useful for solving the full quantum problem. Classical solutions that were useful for the 
problem of operator dimensions were considered in [5, 6] and several other papers. Here we 
consider classical solutions relevant for scattering amplitudes or Wilson loops. 
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A scattering amplitude at strong coupling corresponds to a surface that ends on the 
AdS boundary on a very peculiar polygonal contour [7]. When we consider a color ordered 
amplitude involving n particles with null momenta ki, • ■ ■ , k n we get the following contour. 
The contour is specified by its ordered vertices xi, ■ ■ • , x n , with — x^_ t = fcf , see figure 1. 
The problem becomes identical to the problem of computing a Wilson loop with this contour. 
In fact, we have a "dual conformal symmetry" which acts as the ordinary conformal symmetry 
on the positions Xj [8]. The amplitude has a divergent part and a finite part. The divergent 
part has a structure that is well understood [9]. In addition, a piece of the finite part is 
also known [9, 10]. There is an interesting finite piece which has not yet been computed 
in general. Two loop perturbative computations of this piece include [11, 12, 13, 14], and 
several subsequent papers. The interesting part of the amplitude is a function of conformal 
cross ratios of the Xj. If we have n points we have 3(n — 5) independent cross ratios. At 
strong 't Hooft coupling we can compute this in terms of the area of the minimal surface 
that ends on the polygonal contour [7]. 

Our method will use integrability of the sigma model in the following way. First we define 
a family of flat connections with a spectral parameter 9. Sections of this flat connection can 
be used to define solutions which depend on the spectral parameter 9. With these, we define 
a set cross ratios Y k {9). We find a functional Y system that constrains the 9 dependence of 
the functions Y k . This system has 3(n — 5) "integration constants" which come in when we 
specify the boundary conditions for 9 — > ±oo. We can restate these functional equations in 
terms of integral equations, where the 3(n — 5) parameters appear explicitly. These integral 
equations have a TBA form. Schematically they are 



where the m k and c k are the 3(n — 5) parameters we mentioned above and K rs are some 
kernels. 1 Moreover, the area has an expression in terms of the TBA free energy of the 
system. 



Evaluating Y k at 9 = we get the physical values of the cross ratios. However, we can view 
other values of 9 as a one parameter family of cross ratios which give the same value for the 
area. Thus, changing 9 generates a symmetry of the problem. 

The case involving a six sided polygon was treated in [15] and the octagon, in a particular 
kinematic subspace, was considered in [16, 17]. Using this method, the area is computed 
without finding the explicit shape of the minimal surface. 

Our paper is organized as follows. In section two we recall the connection between sigma 
models which obey the Virasoro constraints and Hitchin equations. In section three we 
discuss the case where the minimal surfaces are embedded in AdS%. This is a warm up 
problem, which is simpler than the general problem. In section four we solve the full AdS^ 
problem. We derive the Y system, the integral equations and the area. We perform some 
checks. We also compute the exact answer for a one parameter family of regular polygons. 
Finally, we present some conclusions. We also have several appendices with useful details. 

1 There are 2n — 10 complex "masses" and n — 5 "chemical potentials" Ck with a precisely reality. 



lognw 



m k cosh 9 + c k + K Ks * log(l + Y s ) 



Area 
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Figure 1: The polygon is specified at the AdS boundary by the positions of the cusps Xj. 
These positions are related to an ordered sequence of momenta kj by kj = Xj — Xj_i. The two 
dimensional minimal surface streches in the AdS bulk and ends on the polygonal contour at 
the boundary. 

2 The classical sigma model and Hitchin equations 

The classical AdS$ sigma model is integrable. This can be shown by exhibiting a one pa- 
rameter family of flat connections. For our problem, it will be convenient to choose this 
one parameter family in a special way which will simplify its asymptotic behavior on the 
worldsheet. In fact, to make this choice we will make use of the Virasoro constraints of the 
theory. This has been explained in detail in previous papers [18, 19, 20, 21, 22, 24]. Instead 
of repeating the whole discussion, we will present a slightly more abstract and algebraic 
version here. 

2.1 General integrable theories and Hitchin equations 

Let us assume that we have a coset space G/H. Let us assume that the Lie algebra Q 
has a Z2 symmetry that ensures integrability. In other words, imagine that the Lie algebra 
has the decomposition Q = TL + K. so that % is left invariant under the action of the Z2 
generator while elements in K are sent to minus themselves. We then write the G invariant 
currents J = g~ x dg. This is a flat current dJ + J A J = 0. We can decompose J in terms its 
components along H and K. as 

J = g- 1 dg = H + K (1) 

When we gauge the sigma model we add a gauge field along H, and we can do local H 
gauge transformations. The equations of motion of the system can be written in terms of 
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the if-gauge invariant currents k = gKg' 1 as d * k = 0. Notice that k are the Noether 
currents of the problem. These equations of motion together with the flatness condition for 
J lead to 

D Z K- Z = = D- Z K Z , [D z , D s ] + [K z , K s ] = (2) 

where D Z X = d z X + [H z , X] (3) 

We can view these as equations for the connection. Once we solve these, we can find a coset 
representative by solving the flatness condition 

(d+ J)^ 1 = (d + H + Kjg- 1 = (4) 

More precisely, we start with a set of independent vector solutions to the equation (d+J)ip = 
0, or t honor malize them, and assemble them into g~ x . These vector solutions are called flat 
sections. The global G-symmetry acts by left multiplication of g and the equations (2) are 
G invariant. The equations (2) are identical to Hitchin's equations after the identification 
$ z = K z , $2 = K z , A = H 2 . These equations are equivalent to the flatness of the one 
parameter family of connections 

K dz 

d + A(Q, with A(C) = -^T + H + ( 2 K- Z dz (5) 

Flat sections of this connection, at £ = 1, give back the group element g~ l . This connection 
differs from the connection that is often written (e.g. in [5]) by a gauge transformation by 
the group element g. Though we will not need it here, let us quote the more usual form of 
the flat connection 

d + a , with a = gAg^ 1 — dgg^ 1 = k z dz{— — 1) + k z dz(( 2 — 1) 

c 2 

We did not find this form of the flat connection particularly useful for our purposes. 

The equations (2) imply that T(z) = Tr[K%] is holomorphic. This is the usual holomor- 
phicity of the stress tensor. For the SO(n + 1)/ ' SO{n) type cosets that we are interested in, 
higher traces of K z vanish, so we do not obtain any other interesting holomorphic quantities. 
In particular, if we are considering a theory obeying the Virasoro constraints, T = 0, then 
we do not appear to get any interesting holomorphic quantities in this fashion. 

2.2 Integrable theories with Virasoro constraints and Hitchin's 
equations 

As we mentioned above, in the case that T = 0, we must work a bit harder in order to obtain 
interesting holomorphic quantities. In fact, it is possible to choose a slightly different form 

2 Actually, to be a bit more precise, we have a Z 2 projection of the Hitchin problem based on G by the 
Z2 symmetry we considered above. Namely, we project on to $ = — s($), and A = s(A) where s is the Z2 
transformation that multiplies the elements of K. by minus one. 
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(or different gauge) for the connection so that we obtain a more interesting Hitchin system. 
This is a small variant of the Pohlmeyer type reduction. For the case with non-zero stress 
tensor this was described in [18, 21, 22]. In the particular case that is going to be of interest 
to us, which is the SO(2,4)/SO(l,5) or AdS$ sigma model this was done in [19, 24, 15]. Since 
we do not want to repeat those derivations here, let us give a more abstract perspective on 
it. 

We consider cosets of the form G/H = SO(n+ l)/SO(n), or SO(2,n)/SO(l,n). Proba- 
bly, similar considerations are true for other cosets but we have not checked the details. We 
now have the Virasoro constraints Trfif^] = and TrfiTf] = 0. We assume that Trfi^i^] 
is generically non-zero. This quantity is the action density or the area element, so it will be 
non-zero for our solutions. We can then think of K z and K s as spanning a two dimensional 
subspace of /C. We consider a generator g in H such that K z has charge +1 and K z has 
charge -1 under q. In other words, we view K z and K z as two lightcone directions in the 
Lie algebra, and q is the "boost" generator. We can further split the lie algebra H according 
to the charges under q. In our case, we have H^°\ and "H*" 1 ), where the superscript 
indicates the charge under q. We then take (5) and make a global gauge transformation by 
( q . We obtain 3 

A = ( q A(- q = \{Kz + tt { ~ x) )dz + H + ((K s + H [ i ] )dz = 5^ + A + ( <$>- z dz (6) 

This is the final form of the flat connection that we will use. We saw that it is a simple 
transformation of the previous one. Moreover, when £ = 1 the gauge transformation is trivial 
and the flat sections of this connection are still giving us the solution g -1 , as in (4). One nice 
aspect is that now P(z) = ^Tr[$^] is a non- vanishing holomorphic current. One can wonder 
why we have a spin four holomorphic current. In general, the integrable theory has higher 
spin conserved currents. These higher spin currents are usually not holomorphic. When the 
stress tensor vanishes, the spin four current becomes holomorphic. In terms of embedding 
coordinates with X 2 = —1, we have P oc d 2 X ■ d 2 X. We also have that Tr[$^] = for r < 4. 
Finally, note that if we start from a general SO(n) Hitchin equation, we can specialize into 
(6) by performing a projection generated by the product of the Z2 transformation we had 
above times a conjugation by (i) q , where q is the U(l) generator we discussed above. This 
combined generator, let us call it r, should then give r($ 2 ) = — i<& z , r($ 2 ) = i$> z , r(A) = A. 
In section 4.1 we give a more explicit form for this generator. 

In our case, we will further use the relation between 5*0(2,4) and SU(2,2) in order 
to write an SU(2, 2) flat connection. If we denote by ip the flat sections of the SU(2, 2) 
connection, then anti-symmetric products of two different sections ip and if)' will give a flat 
section in the vector representation of 5*0(2, 4). Schematically q A = (T A ) a ^if)[ a if)L, where ip 
and ip' are two solutions of the problem in the spinor (or fundamental of S77(2, 2)) and q A 
is a solution in the vector representation of 50(2,4). 

3 In this derivation we have used that tl\ ^ = = H z . This follows from (2) plus the condition that 
K z is non- vanishing (and the g-charges of H and K). 
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Note that the action for the problem, which is equal to the area, is given by 



A = J d 2 zTr S o[2A)[K z K 2 } = 2 J d 2 zTr su{2>2) [K z K,} 

A = J d 2 zTr{$ z ® z ] + total derivative (7) 

The total derivative term is a constant proportional to the degree of the polynomial P 
(and independent of the kinematics). In order to show the last equality in (7) we can take 
the trace of the generator q times the second equation in (2) and use the Jacobi identity. 
(Alternatively, one can show it via an explicit parameterization as in [15].) 

Once we compute this geometric area we can compute the amplitude, or the Wilson loop 
expectation value as 



,2 



Amplitude ~ (W) ~ e ^ a < A = e ^ A 

Here A is the geometrical area of the surface in units where the radius of AdS has been set 
to one. This area is infinite, but it can be regularized in a well understood fashion. The 
central object of this paper is certain regularized area, defined by 



A 



j d 2 z(Tr[<$> z <5>,]-A{PP) l l A ) (8) 



namely, we subtract the behavior of Tr[$ z $ s ] far away. Since (8) is invariant under conformal 
transformations, it is a function of the cross-ratios. When using a physical regulator, the 
area will have additional terms. These additional terms are well understood and described 
in appendix G. 



2.3 Flat sections, Stokes sectors and cross ratios 

In this subsection we recall some facts, which were discussed in more detail in [15]. For the 
amplitude problem the worldsheet is the whole complex plane and P(z) is a polynomial. We 
then study the problem (d+A(())ip = 0. As we go to large z some flat sections ip will diverge 
and some will go to zero. The fact that some diverge means that the worldsheet goes to the 
boundary of AdS space. For large z, the boundary conditions are such that we can simulta- 
neously diagonalize $(z) ~ P(*) 1/4 diag(l,-i,-l,i) and $ s ~ P 1 / 4 (z)diag(l, i, -1, -i). The 
particular relation between eigenvalues is determined by the Z4 symmetry of the problem. 
This determines the large z asymptotics of the four solutions 1 

Va~exp|-r a ^ Pl/ ^ )dZ -i a C J P 1/4 (z)dz\ , a = 0,1, 2, 3 



^A{z) in (6) decays as 1/z for large z and therefore can be dropped when considering the leading asymp- 
totics that determines the Stokes sectors [15]. We will have to keep A when we approximate the cross 
ratios. 
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The problem displays the Stokes phenomenon at large z. This means that the previous 
behavior of solutions is only valid within a given Stokes sector. The number of Stokes sectors 
is determined by the degree of the polynomial. Namely, for large z we have f z dz'P 
z n / A + • • • , for a polynomial of degree n — 4. In order to characterize the problem it is 
convenient to choose the smallest solution in each of the Stokes sectors. This smallest 
solution is well defined up to an overall rescaling. Given any four flat sections, it is possible 
to construct a gauge invariant inner product (ipi, ip 2 , 4>3, V^) = e °' l3 ' y6 i J iai J 2ai J 3'yi J 45- This 
inner product is independent of the position where we compute it. A full solution of the 
problem is given by choosing four arbitrary flat sections tpi, ■ ■ -ip^ where the subindex runs 
over the four solutions, but each of them is a four component spinor. We will use greek 
letters for spinor components and latin letters for labeling different solutions. The target 
space conformal group SU (2, 2) acts on the latin indices, but not on the greek indices where 
the flat connection acts. The spacetime embedding coordinates X 1 are given in terms of 
these solutions at £ = 1. More explicitly, X T T^ b = M^t/w?/^ where M is a fixed matrix 
and a, b are spacetime indices. As we go to large z, some of the solutions diverge. The 
particular combination of solutions that form the two solutions that diverge most rapidly, 
determines a ray X 1 . This maps to a point on the boundary of AdS$ space. We can find this 
point in a convenient way by picking the two smallest solutions which will be s$ and Sj+i, 
if we are between Stokes sectors i and i + 1. Then the spacetime direction is obtained by 
taking 

Xl b cx(il; a ,^,s i ,Si +1 ) (9) 

This determines the direction in which X l ab is diverging. Recall that we can think of the 
boundary of AdS as a projective space, given by six coordinates X 1 , with X 2 = and 
X 1 ~ XX 1 . Thus the diverging solution determines a point X in projective space, which is 
the same as saying that it determines a point on the boundary of AdS. 

The index % labels the cusp number. We can form quantities of the form X 1 ■ X^ oc 
(sj, Sj + i, Sj, Sj+i). Finally, cross ratios are given by quantities of the form 

Y = X { -X j X k -X l = (s i ,Si+i,s j ,s j+1 )(s k ,s k+1 ,si,si +1 ) 
vki x i -X k X l -Xi (s i ,s i+1 ,Sk,s k+1 )(si,si +1 ,s j ,s j+1 ) 

The cross ratios do not depend on the normalization of each of the Sj. 

These cross ratios are functions of the spectral parameter (. In what follows, we will 
choose a convenient basis of cross ratios and study their ( dependence. We will write an 
integral equation determining the values of the cross ratios as a function of (. Finally we 
will express the area in terms of certain integrals of the cross ratios over (. 



3 Minimal surfaces in AdSs 

We first consider minimal surfaces that can be embedded in an AdS% subspace of AdS§. 
This is a simpler problem that illustrates the method that we will use in the AdS^ case. The 
reader that is only interested in the AdS$ case can jump directly to the next section. 
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3.1 AdSs preliminaries 



When the surface can be embedded in AdS^ the problem simplifies and it reduces to a Z 2 
projection of an SU(2) Hitchin problem. The derivation of this fact is rather similar to what 
we discussed above and was treated in detail in [17]. We will not repeat the derivation, but 
we will state the final results. We have a polynomial p = |Tr[$^] whose degree determines 
the number of cusps. 5 We now have SU(2) quantities § z , A, $ f which are in the adjoint 
of SU(2) and we have the Z 2 projection condition $ 2 = —U§ z U~ l , $ 2 = —U&gU^ 1 and 
A = UAU~ l where U = 03 is the usual Pauli matrix. This restricts the components of A and 
$ that are non-zero. General SU(2) Hitchin problems were studied in [25] and we are now 
considering a special case of their discussion, though we will rederive some of their formulas 
in a different way. We study sections of the flat connection which obey 

(d+^ + A + <S>- z dz CWO=0 

The Z 2 symmetry relates solutions V'(C) with different values of the spectral parameter. 
Namely, if ip(C) is a flat section with spectral parameter (, then 77(C) = Utp(e l7r Q is a 
solution of the problem with spectral parameter (. We can track how the small solutions 
change as we change ( by looking at them in the large z region. In a given Stokes sector 

the small solution contains a factor behaving as e « ~ e c , where n is determined 
by the degree of the polynomial and is equal to the number of cusps of the polygon, n is 
even. There are n/2 Stokes sectors and thus n/2 small solutions Sj. As we change the phase 
of (, the ray in the z plane where this solution is smallest rotates accordingly. In particular, 
if we start with the solution which is the small in the ith Stokes sector, we find that 

Si(e 2m () s i+2 (C) an d Usi(e l7T () oc s i+ i((). Note that the solutions do not come back to 
themselves after a shift by e 2m (. We can choose a solution s\ in the first Stokes sector and 
define all others as Sj = U^^Si^e^Q. Then, as we go around we have that sa +1 = A(Qsx. 
A(() can be set to one when n/2 is odd. When n/2 is even it has a simple form that we will 
discuss later. 

The full connection with spectral parameter is an SL(2) connection and thus we can form 
an SL(2) invariant product (ipip') with two solutions. Now we have that 

( Si , Sj )(eH) = (si + i,s j+1 )(C) (11) 

We can normalize S\ so that (si,^) = 1. Then (11) also implies that (sj,Sj + i) = 1. 
We can form cross ratios by forming quantities like 

x.MO = ! S "^M (12) 

These quantities do not depend on the arbitrary normalization of the Sj. By construction 
they are also invariant under the conformal symmetries of Ad S3. They can be related to 
the conformal invariant cross ratios formed from the positions of the cusps of the polygon. 



^Note that P oc p 2 , where P — jTr($^) i n the AdSs polynomial discussed in the previous section. 
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Figure 2: Spacetime positions of the cusps for a polygon that is embedded in R > , which is 
the boundary of AdS^. The positions of the cusps are given by a set of n/2 values xf and a 
set of n/2 values x~[ . 



Recall that a polygon in AdS^ is given by n/2 positions x\ and n/2 positions x~ , see figure 
2. We can form spacetime cross ratios from the positions of the points xf. These spacetime 
cross ratios can be expressed in terms of the cross ratios in (12) as 

xmiC = i) = -Vt (13) 
xm(C = i) = ^ (14) 



3.2 The functional Y-system 

We will now derive a set of functional equations for the inner products, or Wronskians, 
(si,Sj)(() made out of two small solutions of the linear problem. The starting point is the 
Schouten identity, (sj, Sj){s k , si) + (s i; si){sj,Sk) + (s,, Sk)(si, Sj) = 0, applied to a particular 
choice of small solutions: 

(Sjfc+i, S_ fc )(s fc , = (Sfc+i, S_ fc _i)(Sjfc, S_ fc ) + (Sjfe, Sfc + i)(s_jfe_ 1 , s_ fc ) . (15) 

In our normalization the last two brackets are equal to one. Using (11) we see that this 
identity becomes the SU(2) Hirota equation 

r s +T s - = r s+1 r s _ 1 + i, (16) 

where : T 2k +i = (s-fc-i, s k+1 ) , T 2k = (s_ fc _i,s fc ) + (17) 

or more uniformly T s = (sq, s s+ i)(e _ ^ s+1 - )7r / 2 C). The superscripts ± indicate a shift in spectral 
parameter, f ± = f{e ±% ^C/). Actually, from (15) we get (16) for s = 2k. For s odd we need to 
start from a slightly different choice of indices in (15). T s is non-zero for s = 0, . . . , n/2 — 2. 
Finally, we introduce the F-functions Y s = T s _iT s+1 . Being a product of two next-to-nearest- 
neighbor T-functions, the F-functions are non-zero in a slightly smaller lattice parametrized 
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bys = l,...,n/2 — 3. The number of Y- functions coincides with the number of independent 
cross ratios. 

The Hirota equation (16) implies the Y-system for these new quantities: 

Y s + Y- = (l + Y s+1 )(l + Y s ^). (18) 

These equations are of course not enough to fix the Y-functions. After all they came from a 
trivial determinant identity without any information on the dynamics! To render them more 
restrictive we need to supplement them with the analytic properties of the Y-functions. 
This will then pick the appropriate solutions to these equations. Furthermore, to make 
these solutions useful we must relate them to the actual expression for the area. Before 
considering these points let us comment on some general properties of Hirota equations and 
their corresponding Y-systems. 



3.3 Hirota equation, gauge invariance and normalization of small 
solutions 

The general form of the Hirota system of equations - which generalizes the £77(2) case 
derived above - is a set of functional equations for functions T ajS (£). 6 The indices a, s take 
integers values and can be thought of as parametrizing a two dimensional lattice. At each 
point of this lattice we have a function T a s (Q of the spectral parameter (. Then, for each 
site o = (a, s) we have an Hirota equation 

T+T- = T+.T.+ + T t T 4 (19) 

involving the function at that site and the four T-functions at the four nearest-neighbor sites, 
= T 0jg+ i, = T a+ i :S , etc. Recall that = T (e ±t7T ^ 2 Q. This equation has a huge gauge 
redundancy 

a,/3=± 

where g a p{() are f° ur arbitrary functions. It is therefore instructive to construct a set of 
gauge invariant quantities 

Y = 7 ^ or Y a , s = Z a ' s - l7 }' s+1 (20) 

It is instructive to think of the gauge invariant quantity Y Q as a field strength made of the 
gauge dependent gauge field T Q . Suppose the T-functions are non-zero in some rectangular 
domain in the (a, s) lattice. 7 At the edges of the rectangle either the first or the second term 



6 Typically these relations arise in the study of quantum integrable models and describe the fusion relations 
for the eigenvalues T a s of transfer matrices in rectangular representations parametrized by Young tableaus 
with a rows and s columns [23]. 

7 Strictly speaking the T-functions can not be non-zero only inside a finite rectangle: by analyzing Hirota 
at the upper right corner (a*, s*) of the rectangle we would conclude that T a *,s* = 0. This would then imply 
that the neighbors of this point, T *_ liS « = T a ». s »_i = which will then imply that T a »_2.s* = T a * lS *-2 = 
etc; at the end we would be left with T QiS = everywhere. What we can have, for example, is T a ^ s ^ in 
a rectangle and on the two infinite lines containing the upper and lower edges of the rectangle, see figure 3. 
At these lines T ajS are trivial (pure gauge) but they are non-zero. This is what we mean in the text. 
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in the right hand side of (19) is zero. We are left with a discrete Laplace equation for (the 
logarithm of) the T-functions and therefore they become pure gauge. At these boundary 
points the Y-functions are trivial (either zero or infinity) as expected from the analogy. The 
Y-functions are non-trivial in a smaller rectangle obtained from by removing the first and 
last columns and rows of the original domain. 

The Hirota equation 19 then translates into the Y-system 

y + y- = (i + y^)(i + y_) 

y t n (i+F t )(i+n) 

for these gauge invariant quantities. Different domains in a, s where Y's are nontrivial to- 
gether with different boundary condition and analytic properties describe different integrable 
models. 

In the treatment of the previous section we considered the case where the T-functions live 
in a finite strip with three rows and n/2 — 1 columns, where n is the number of gluons, see 
figure 3. The functions denoted by T s in that section are the T-functions in the non-trivial 
middle row, T s = T\ iS . Similarly Y s = Yi jS . The T-functions introduced in that section are 
inner products of small solutions and are therefore sensitive to their normalization. This 
arbitrariness is a manifestation of the gauge freedom in Hirota equation. The normalization 
(si, s i+i) — 1 corresponds to the gauge choice where 

To,2k = ( s -fc-l; s -k) ) T),2Ai+l = ( s -k-2, s -k-l) + > 

^~2,2fc = (Sk, Sjfc+l) > T 2: 2k+1 = (Sk, Sk+l) + , 

are gauge fixed to one. We could of course opt not to fix a normalization for the T-functions 
but then we should use the gauge invariant combination (20) when defining the Y-functions: 

\r T.,2fc-lT.,2fc+l _ \S-k, Sk)\ s -k-l, s k+l) / 00 \ 

Y 2k — — 7p — 7p — 7 T7 r 

I 0,2k I 2,2k W-fc-l, S_k/\Sk, Sk+l) 



Y 



2k+l 



Tl > 2fcTi i 2A:+2 
T),2fc+lT2 ) 



2k+l 



(S-k-1, Sfc)(-S-fc-2, -Sfc+l) 
_(s_ fc _ 2 , S- k -l)(Sk, Sk+l) 



We see that they are now manifestly independent of the choice of the normalization of the 
small solutions. At spectral parameter £ = 1, or £ = e ±t7T / 2 they yield physical space-time 
cross-ratios as in (13) (14). 

A particularly interesting quantity is 

B(Q = Yn / 2 ~ 3 Yn/2-l Y n /2-H 
Y„/2-5 Y n /2_9 Y n /2_13 

It follows from the Y-system equations (18) that B + B~ = 1. B(Q is constructed from the 
Y-functions and is therefore gauge invariant. Using the definition of the Y-functions we 
see that B(() is given by a bunch of boundary T-functions. In our normalization all these 
functions except for the rightmost one are gauged to one. We find therefore B = T\ >n /2-2 = 
(si,s n /2) (e~ t7T ( n / 2+1 ^ 2 () . This means that B(Q is the function that governs the monodormy 
Sn/2 = -B(Ce i7r{n/2+1)/2 )s of the small solutions. 
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Figure 3: In this figure we have summarized the structure of the T's and the V's in a gauge 
where we simplified the T's that can be simplified. The small solid black dots represent 
non-zero T-functions. They are equal to one unless a = 1 and s — 1, . . . , n/2 — 2. At the 
rightmost point in this line we have T l ra / 2 _ 2 = B where B is a function which cannot be set 
to one only in the case that n/2 is even. In fact, in our case it is B = —e m ^ +m< ^. This is the 
function that governs the monodormy s n / 2 = — B(^e ln( - n ^ 2+1 ^ 2 )sQ. The ^-functions are finite 
in the points indicated by fat gray shaded balls. At all other points they are either zero or 
infinity. 



3.4 Analytic properties of the y-functions 

For finite values of ( it is clear from (17) that the T s are analytic functions of (, for ( ^ 
0, oo. Generically, they will not be periodic under ( — > e 2m (. In general, the Y's will be 
meromorphic functions. However, in our case, since we can choose to set the denominators to 
one, we see that the Ys have no poles and are thus analytic away from £ = 0, oo. For £ — > 
and ( — > oo they will have essential singularities. In this section we analyze the behavior in 
these two regions. 

When ( — > we can solve the equations for the flat sections by making a WKB ap- 
proximation, where ( plays the role of h. This is explained in great detail in [25], here we 
will summarize that discussion and apply it to our case. The final result is that, for an 
appropriate choice of the polynomial p, we have the standard boundary conditions in TBA 
equations. We will later discuss what happens for more general polynomials. 

We are considering the equation 

(d + + A + C ®zdz^J s = 

When ( — > 0, it is convenient to make a similarity transformation that diagonalizes $ 2 — > 
v /pdiag(l, —1). The solutions in this approximation go like exp J y/pdz^j times constant 
vectors. The WKB is a good approximation if we are following the solution along a line of 
steepest descent. This is a line where the variation of the exponent, is real, lm{^p{z)z/Q = 
0. This condition is an equation which determines the WKB lines. Through each point in 
the z plane we have one such line going through. At the single zeros of p we have three lines 
coming in. The WKB approximation fails at the zeros of p (which are the turning points). 
From each Stokes sector we have WKB lines that emanate from it. These lines can end in 
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Figure 4: Sketch of the pattern of WKB lines for various phases of (. The crosses denote 
the various zeros of p(z). The numbers indicate the various Stokes sectors. The black thin 
lines end on the zeros and separate different classes of WKB lines. The thick colored lines 
are the WKB lines that we use to evaluate cross ratios. Here we have indicated only the 
ones used to evaluate Y 2 and Y 3 . Finally, note that by setting the phase of ( to e ln ^ we have 
WKB lines that enable us to evaluate all the Y s . 

other Stokes sectors or, for very special lines, on the zeros of p. If a line connects two Stokes 
sectors, say i and j, then we can use it to approximate reliably the inner product (si,Sj). 
This estimate is good in a sector of width tt in the phase of (, centered on the value of ( 
where the line exists. As we change the phase of ( the pattern of flow lines changes. It also 
changes when we change the polynomial p. We first select a polynomial with all zeros along 
the real axis and such that p(z) > for large enough values of z along the real line. 

With this choice the pattern flows for WKB lines is shown in figure 4, for some values of 
the phase of (. The WKB lines ending on zeros separate regions where lines flow between 
different Stokes sectors. In our problem we have some inner products evaluated at ( and 
some and i(. The flows for ( = 1 and ( = i are displayed on the top of figure 4, and they can 
be used to evaluate the various inner products. Alternatively, we can set ( = e l7r//4 , evaluate 
them all and then continue them from this region. The resulting flow pattern is sketched on 
the bottom of figure 4. 

Using those flow patterns it is a simple matter to evaluate various inner products. It 
turns out that the inner products in the definitions of the Y-functions (22) combine to give 
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Figure 5: Cycles along which we need to integrate yfpdz in order to determine the asymptotic 
form of Y s . By being careful about the sheet selected by the various small solutions one can 
determine the cycle orientations shown here. 

a contour integral around a certain cycles. See figure 5. Thus, each Y s , is estimated by the 
integral of y/p along a cycle j s . We can call 




and the corresponding Y functions have the small ( behavior 

log Y 2k — H , log Y 2k +i 7— H 

c < 

In figure 5 we display the cycles corresponding to each of the Y s . It is convenient to define 
the parameters m s via 

m 2k = -2Z 2k m 2 k-i = -2Z 2k+1 /i (24) 

For our choice of polynomial the m s are all real and positive. In order to check the positivity 
of the m s we need to be careful with the choice of branch when we evaluate the cross ratios. 
The two branches correspond to the two eigenvalues of $ z and differ by an overall sign. 
Taking the same cycles but on different branches is equivalent to changing the sign of Z. 
The correct branch is determined by the behavior of the various small solutions, each of 
which goes like e ± ^ s/v dz . After taking this into account we can check that the m s are indeed 
positive for the polynomial we chose. 

A similar computation at large ( gives a similar result, with logY^fc ~ (Z 2k . Thus, we 
have shown that all the Y functions have the asymptotic behavior 

log Y s = —m s cosh 8 + ■ • • 

for large 8, ( = e e . 

Furthermore, this behavior is good over a range of (— 7r, 7r) in the imaginary part of 8 s . 
The reason is the following. For each Y s , the region of Im{6) = corresponds to the center 
of the region where the WKB line exists. Furthermore, the corresponding WKB lines exist 
for a sector of angular size 7r around this line. In addition, we have mentioned that the WKB 
approximation continues to be good for a further sector of tt/2 on each side. In fact, this is 
more than enough for deriving the integral equations. 

8 Recall that the functions are not periodic in lm(8). 
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3.5 Integral form of the equations 



The analytic properties derived above together with the functional equations (18) uniquely 
determine the Y-functions. However, for practical purposes - specially for numerics - it is 
useful to have an equivalent formulation of these Y-system equations in terms of TBA-like 
integral equations. 

To derive them we follow the usual procedure which we briefly review for completeness. 
We note that l s = log(Y s /e _ms cosh6 ') is analytic in the strip |lm(0)| < n/2, vanishes as 9 
approaches infinity in this strip, and obeys 



e + /; = iog(i + r s+1 )(i + K 



8-1 j 



which is nothing but the logarithm of the Y-system equations. Now we convolute this 
equation with the kernel 



K{9) ' 



2n cosh 9 



to get 

K * (it + r s 



dy l 8 (y + iir/2) + l s (y -in/2) _ £ dy l s (y) _ ^ 



2tx cosh(x — y) J 2iri sinh(x — y) 



where 7 is the rectangle made out of the boundaries of the physical strip together with two 
vertical segments at Re(#) — > ±00. In order to be able to add these extra segments to the 
integral thus making it a contour integral, it was important to use the l s instead of log Y s . 
This is why this quantity was introduced in the first place. Furthermore, in the last step we 
used the fact that l s has no singularities inside the physical strip, this is an important input 
on the analytic properties of the Y-functions. Rewriting l s in terms of Y s leads therefore to 
the desired form of the integral equations: 

logr s = -m s cos\i9 + Ki<\og(l + Y 3+1 )(1 + Y S _ 1 ) . (25) 

For a given choice of masses m s , the solution to these integral equations is unique and a 
basis of cross ratios can be read from evaluating the Y s (9) , s at 9 = 0. These equations are 
of the form of those appearing in the study of the TBA [1] for quantum integrable models 
in finite volume. Furthermore, as explained in the next section, the (regularized) area of the 
minimal surfaces turns out to be given in terms of the Y-functions as the free energy of the 
corresponding integrable model. 

Up to now we have discussed the case where the zeros of the polynomial are along the 
real axis. Let us briefly discuss what happens as we start moving the zeros of the polynomial 
away from the real axis. Notice that the functional Y system equations (18) do not depend 
on the polynomial. Thus, these equations continue to be true, regardless of its form. What 
changes are the asymptotic boundary conditions. 

Let us consider first the case were we start from the above polynomial and we move the 
zeros around a little bit. Then the above derivation of the asymptotic form of the Y functions 
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goes through with only one change. Namely the quantities Z s and m s are now more general 
complex numbers and the asymptotic behavior is 

logF, ~ , for c -»■ o 

logF s ~ "'Y^ ' for C^oo 

In this case, when we derive the integral equation, it is convenient to shift the line where 
the Y functions are integrated to be along the direction where m s /( = \m s \e l(ps /£ is real and 
positive, which also makes fh s ( real and positive. Then, defining Y s (6) = Y s {9 + iip s ) we find 
that the integral equations have the form 

logy = -\m s \ coshtf + K s , s+1 *log(l + y+i) + tf.,._i*log(l + y s _x) (26) 

where now K ss i = 1/ cosh(# — & + i<f s — i<f s ')- 

As long as \(p s — <^ s +i| < 7r/2, the integral equations conserve the form that we have 
derived. If we deform the phases beyond that regime we will have to change the form of 
the integral equations by picking the appropriate pole contributions from the kernels (which 
become singular for \ip s — (p s+ i\ = tt/2, 3tt/2, . . .). Of course, the integral equation changes 
but the y's and therefore the area are continuous. The pattern under which the integral 
equation changes is explained in appendix B. 9 This is the wall crossing phenomenon d iscussed 
in [26, 25]. 

These integral equations are a special case of the general case discussed in [26]. In fact, 
the equations in [26] are true for an arbitrary M = 2 theory, and a Hitchin problem is just 
a special case. Due to the Z 2 projection we have that the X 1 quantities in [26] obey the 
additional property AL 7 (£) = C). Using this, we can easily map the kernel in [26] to 
the -K- found here. 



3.6 Area and free energy 

As we mentioned above, the interesting part of the area is given by the integral 

A = 2 J d 2 zTr{$ z $ 2 } (27) 

By definition, the area is independent of 

It is convenient to think again about the small ( regime and the WKB approximation 
that we did for small (. In fact, we can improve on the WKB approximation and find the 
next couple of terms by systematically expanding the expressions for the inner products. We 
take complex masses but with small enough phases so that the WKB approximations that 

9 A very similar kind of manipulations is outlined in section 4.7 when we explain in detail how to compute 
the F-functions in the AdS$ case for large values of the imaginary part of 9. 
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we did before continue to be valid, with the same cycles. 10 . The final At this point it is 
convenient to use slightly different functions defined by 



in order to undo the shifts in (22). With these definitions we find that 

"f A 



logY 



c 

A = xdz 



a + ( (t u 4 



(28) 



(29) 



x 



p(z) 



where u is an exact one form. Here a given by the diagonal components of the connection 
A. In our case a = due to the Z 2 projection. (But even if a were non-zero, it would not 
affect what we say below.) We know that u is exact because we can deform the contour and 
logK, should not change. It has a u z and a u s component. For our purposes, it will only be 
important to compute the w 5 component which is 



(30) 



where the index i is not summed over. In other words, we get the diagonal components of $2 
and we get the first or second diagonal component depending on whether we are on the first 
or second sheet of the Riemann surface. That is, u is a one form on the Riemann surface, 
not on the z plane. In the basis where <3> 2 is diagonal, we can thus rewrite (27), using (30), 
as 



A 



A Am 



UK. 



A ® 7/ 



where Y are a basis of cycles 11 . We will take this basis to be the basis of cycles that gives 
the WKB approximation to the Y s , see figure 5, and w rs is the inverse of the intersection 
form of the cycles. The matrix of cycle intersections can be read off from figure 5 and it is 
summarized in figure 6. 

We can also compute the small ( behavior of Yf, by expanding the integral equations (25). 
We get 



iogY P = zvc + C 



5> 



1 



Til 



c c 



7 iog(i+n: 



(31) 



It turns out that 9 rs is given by the intersection form of the cycles involved. This follows 
from the general theory in [26], but it can be easily checked in this case by examining the 



10 In general, the cross ratios that have a simple WKB approximation will change as we change the phase 
of the masses beyond a certain point, see [25]. For us, it is enough to do the derivation for some range of 
masses. Then we can analytically continue the final formula, as explained in appendix B. 

11 This formula looks suspicious because the left hand side is infinite while the right hand side is finite!. 
Here we have implicitly used a regularization which puts a cutoff in the w plane for large values of \w\, where 
dw = yfpdz. We have then subtracted the same integral but with a polynomial whose zeros are all at the 
origin. This procedure works well when n/2 is odd. 
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Y 1 +Y 2 +Y 3 +Y 4 +Y 5 +Y 6 +Y 1 +Ys+Y 9 +Y l0 -+Y n -+-Y l2 -~Y l3 

Figure 6: Intersection form 9 r,s for all the cycles associated to the Y functions. This is 
computed from the cycles in figure 5. If an arrow points from Y s to Y r , then we have 
(ls,Jr) = I? otherwise the intersection vanishes. 

integral equations (25) and remembering that the Y s differ by simple shifts in the argument 
from the Y s . 

Thus, we obtain that A = A perio ds + Af ree with 

^■periods %W rs Z Z 



Afree 



-;£/ff>«< i+? -> (32) 

Afree = ^ / — \m 8 | cosh 9 log(l + Y s ) (33) 



with Y s (9) = Y S (9 + icp s ) , m s = \m s \e 

In order to obtain (33) we have averaged the result from (32) with the result we obtain 
from the large ( expansion. The fact that large £ and small £ should give the same answer 
translates into the statement that the total momentum of the TBA system should be zero. 
The explicit form of A perio d s in terms of the masses is given in (96). 

This derivation has assumed that n/2 is odd, because we said that 9 rs was invertible. If 
n/2 is even, then we start from n/2 + 1 and we take away one zero of the polynomial. Then 
the result contains two pieces, one piece has the form of Af ree discussed above and the other 
contains an extra term that was discussed in detail in [17]. 

Also, in this derivation, we have assumed that the intersection form of the cycles asso- 
ciated to the Yy that appear in the integral equation is invertible. While this is true in our 
case, it would cease to be true once we cross walls and we get extra cycles [26]. One can 
slightly modify the above derivation and the final answer continues to be (33), see appendix 
B. 

Note that, in the end, we do not need to know the polynomial, or the Riemann surface, 
or the cycles. That is only needed for the derivation. Ultimately, everything is expressed in 
terms of the (complex) masses m s appearing in the integral equations, see (33), (96). 



3.7 The octagon, or n = 8 

Here we rederive some of the results in [17] from this point of view. In this case there is only 
one Y function and the functional equation is Y + Y~ = 1, whose solution is just Y = e z ^ +zc -. 
The free energy is 

A free = ^- [ ^2|Z|cosh#log(l + e- 2|z|coshe ) 
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which agrees with what was called A sin h in [17]. 

The full result in [17] contains an extra piece A extra which is related to an extra com- 
plication that arises in the case that n/2 is even. In this case we will also need the Hirota 
variable T\. 



4 Minimal surfaces in AdS§ 
4.1 AdSz preliminaries 

As we mentioned in section two, the worldsheet theory describing strings in AdS$ can be 
reduced to a Z 4 projection of an Si7(4), or SU (2,2), Hitchin system. 

After some gauge choices we can represent the action of the Z4 in the following way 

/ 1 \ 

10 

1 

\ -1 / 



r(X) = -CX'C- 1 , C~ 



(34) 



where A is $ or A. Recall that we will be imposing the projection conditions r($ z ) = — z<3> 2 , 
r(A) = A and r(& 2 ) = 

This Z 4 symmetry relates solutions to the problem at i( with solutions to a related 
problem at (. More precisely, they relate solutions to the "inverse" problem at (. More 
explicitly, we have 

<H0 = c-V«) 

where V>(C) is a flat section of (d — A 1 )^ = 0. Note that the bar does not denote complex 
conjugation. Given a solution ^(C) of the straight problem and a solution of the inverse 
problem, one can form an inner product of the form This inner product can be 

computed at any point on the worldsheet, and it will be independent of the point where it 
is computed. 

Another property that we will use is the following. Imagine we start with three different 
solutions of the linear problem ipi, ip2, ^3- Then the following combination is a solution of 
the inverse problem 

^23 = C^^Vl/Ji ^2/^3/33 =^lA^ 2 A^ 

where the last equality is simply the definition of the wedge product. This gives something 
interesting when it is applied to small solutions Si for three consecutive Stokes sectors 

Sj_i A 5^ A s i+1 ex Si 

In other words, this product of three solutions gives us a solution to the inverse problem 
that is small in Stokes sector i. This follows from the asymptotic behavior of the solutions 
in Stokes sector i. 
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We show in appendix C that we can choose normalizations for all solutions, Sj, s,, so that 
the following equalities are true 

1 = (Si, Sj+l, Si + 2, Si + z) 

Si = Si_i A Si A s i+1 (35) 

s i+1 (c) = c-h^o 

s i+1 (C) = C T s t (e^ 2 () 
Using these formulas, plus identities involving e symbols, it is possible to show that 
(s k ,s k +i,Sj,s j+ i)(0 = (sk-i,Sk,Sj-i,Sj) (e I7r/2 C) 
(s i ,s fc ,s fc+ i,s fc+2 )(C) = (s i ,s j _i,s j _2,s fc )(e m/2 C) (36) 
(%,s fe+1 )(C) = (^_ 1)S ,)(e-/ 2 C) 
Where the last two lines are stating the same result in two alternative notations. 

Finally, if the problem involves n Stokes sectors, we expect that Si +n oc s», where Si +n 
has been obtained by going around the z plane and normalizing the solutions via (35). We 
will not need the proportionality constant to derive the Y system. In fact, one can calculate 
it from the Y system. However, it is possible to compute it just from the behavior of the 
solutions at infinity, see appendix C. The result is that for n odd we can normalize the 
solutions so that s n = Sq. When n is even this is not possible. When n = Ak one has 

-i^+WriC x -7+ — IWnC + — WqC /n>7\ 

s n = /ie« us s , s n -i = -e< s n _ 2 = /ie < s_ 2 (37) 

A 4 

where wq is a constant. For n = Ak + 2, wq = and only \x is allowed. 



4.2 The AdS 5 Y-system. 

The basic identity which we will use to derive the Y and T-system functional equations for 
the spectral parameter dependent products, or Wronskians, is a particular case of Plucker 
relations: Let (x\, . . . ,xn) be the determinant of the N x N matrix whose columns are N- 
vectors Xi. Then, for iV + 2 generic vectors x\, . . . , xn, Hi, Un we have the following identity 

(xi, * , x N ) (yx, * , y N ) = (yi, * , x N ) (x u * , y N ) - (y N , * , x N ) (x u * , y x ) , (38) 

where * = x%, . . . , xjv-i. 

For us, N = 4 and (x\, *, X4) = (si, Sj, Sk, si) are determinants of 4 x 4 matrices composed 
from four sections of the linear problem. Using the small solutions {s m , s m+ i, s m+ 2, s m +3, S-i, «o}, 
{s_ 2 , s m+1 , s_i, s m , s , s m+2 } and {s_ 2 , s_i, s , Si, s m+ i, s m+2 } for {a;i, x 2 , x 3 , x 4 , y u y 4 } we ob- 
tain three Hirota relations 

T'+ rp — rp rp I rp rp 

l,m I 3,m — 1 l,m-l-*3,m+l + J 0,m-L2,m , 

T'+ T^— T -1 1 T T 

2,m 1 2,m ~ 1 2,m+l 1 2,m-l t J l,m J 3,m > 

3,m J l,?Tt — 1,?72+1 3,m— 1 i" 4,m-^ 2,m • 
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In a more compact notation these read 



"^ffl m -^~4— a,m+l -^~a,m— 1 "I - -^~a+l,m T a _ l,m j (^9) 

where a = 1, 2, 3 correspondingly. The T-functions are given by 

^0,m(C) = ( s m, S m+ i, S m+ 2, S m+ 3)' m ^ = 1 , T 4]?n (£) = (s_2, S-l, Sq, Si)' m ^ = 1 , 

and most importantly 

Tl,m(C) = ( s -2, S-l) So, S m+ i}' ™' , 

T 2 , m (C) = (S-i,So,Sm + l,Sm+2} hm " 1] , (40) 

T 3 , m (C) = S m , S m+ i, S m+2 ) [_m] • 

Here the superscripts ± and [m] indicate shifts in the spectral parameter, f ± = f(e ±l ^() 
and = f(e l ^ m Q. Note that these shifts are half the ones appearing in the Ad S3 case. 12 
When deriving the Hirota equations from the Plucker identity we use the relations (36). 

Equation (39) is an exotic form of the Hirota equation (19). It is exotic because of the 
appearance of T 4 _ 0jS instead of the usual T a>s . 

For s — — 1 and s = n — 3 we see that T a s = 0. This follows simply from the defmtion of 
the T-functions together with Sk+ n oc Thus, the functions T a>s live in a strip of five rows 
(a = 0, 1, 2, 3, 4) and n — 3 columns (s = 0, 1, . . . , n — 4), see figure 7. As explained in section 
3.3, the Hirota equation has a huge gauge symmetry. The same is also true for the more 
exotic form considered here. 13 Our choice of normalization (siS i+ iSi + 2Si +3 ) = 1, corresponds 
to the gauge choice where the T-functions at the left, top and bottom bondaries of the strip 
are set to one, 

T a ,o = T , s = T 4 , s = 1 . (41) 

At the right boundary of the strip, in our normalization (35), the T-functions are related to 
the formal monodromy which arise when comparing Sj with s i+n (37), 



(42) 



J- 1,71-4 — ^ |e <• I ,-l2,n-4 - |C <• I ,J3 jn _4 — /ile « 

Next, we introduce the "K-functions 



J a+l,m-' a— l,m 

Being composed from next-to-nearest-neighbor T-functions, the y-functions y ajS are finite 
in a slightly smaller lattice parametrized by a = 1,2,3 and s = 1, . . . , n — 5. The number 



12 Whether we are discussing the AdS^ or the AdSc, case will be clear from the text. We hope that the 
difference between the shifts in the two cases will not cause a confusion. 

13 For example, if g(a, s, 9) is a gauge symmetry transformation of Hirota (19), then g(a, s, 0)g{A — a, s, 9) 
is a gauge symmetry of the more exotic form (39). 
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Figure 7: Strip where the T and Y-functions live in the AdS§ case. Small solid black dots 
represent T-functions. At the boundary the T-functions are equal to one except at the three 
nodes in the right boundary; there they take the values indicated in the figure. In all the 
points of the boundary the Y-functions are either zero or infinity. They are non-trivial in 
the smaller domain indicated by the fat shaded gray circles. 



of y-functions coincides with the number of independent cross ratios. The Hirota equation 
then implies the following Y-system for these quantities: 14 



y- y+ 

2,m 2,m 

yt y + 



y 



2,m 



3,m 



,m+l ,m—l ) 

(i + y lim )(i + y 3 , m ) 

1+^3 ,m+l 

)(i + y 1 , m _ 1 ) 

;i + y ,m+l ,m— 1 ) 



(44) 



y 



2,m 



or in a more compact notation 

(i + y ,m+l )(l+^4 — a,m— 1 



y- y+ 

a,m 4—a,m 

V V 

1 a+ljTn 1 a—l,m 



(1 + io+l,m)(l + ^a-l,r. 

To recover (44) we need to use that y 0)S = Y^ s = oo. 



a = 1, 2, 3 , 

s = 1, . . . , n — 5 



(45) 



4.3 Analytic properties of the y-functions 

To derive the integral form of the Y-system equations it is important to identify the large 
9 asymptotics. They are fixed by a WKB analysis. The method is very similar to the one 
used for the AdS^ case, but a bit more involved. We will leave the details for appendix E 
and state here the final results. 

14 Similar exotic ^-systems recently appeared in a very different context [27] 
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We choose the polynomial P to be such that all zeros are on the real axis and P(z) > 
for sufficiently large z. Then the large 9 behavior of the Y-functions is 

log Y\ yS — > —m s cosh 9 — C s ± D s , 9 — > ±oo 

log Y 3jS -> —m a cosh 9 + C s =F D s , 9 — >■ ±oo (46) 

logY 2jS — — m s cosh 9 , — > ±oo 

where 6 1 = log(£). The constants C s and D s arise from the components of the connection A 
that survive the projection. For loops in signature (1,3) or (3, 1) the D s 's are real while 
the C s 's are purely imaginary, see appendix F. In fact, we have the general reality condition 15 

CMC))* = *WVC*) 

In particular, for large 9 = log (, we see that 

[log(Y ljS /Y 3 , s ) (+oo)}* = -log(F M /y 3 ,s) (-oo) 

which indeed implies what we said above regarding D s and C s . It turns out that the m s can 
be promoted to complex constants by changing the position of the zeros of the polynomial 16 . 
These n — 3 complex constants, together with the purely imaginary C s , constitute the 3(n — 5) 
parameters of the problem. 

4.4 TBA equations in AdS 5 

To derive a set of integral equations from the functional Y-system we follow again the same 
route as in the AdS 3 case. As explained in the Ad S3 case, a big advantage of the integral 
form of the Y-system equations is the straightforward numerical implementation of these 
equations. We first consider the case where all masses are real and positive. We then 
introduce a set of functions l a , s = log(Y ajS ) + m O)S cosh0 which are meromorphic in the strip 
—71 /A < Im(9) < 7r/4 and bounded as we approach infinity inside this strip. Let us first 
assume that these functions are actually holomorphic in the strip and then we will mention 
what happens when there are poles. Equally important they obey 

C + l 4-a,s ~ l a+l,s ~ L~l,s = log Y+ s + log Y^_ a s - log F a+M - log Y a _i, s . 

The right hand side of this equality is the logarithm of the left hand side of the Y-system 
functional equations (44) derived above. Now we go to Fourier space where we have 

where J 7 denotes the Fourier transform. When writing this relation we are making use of 
the analytic properties of l a>s mentioned above. The Y-system equations can then be cast as 

A ! \T(1 \t \ T (\ (1 + ^a,s+l)(l + ^4-a,s-l) \ / x 

A aa ,(u)T{l a , tS )(u) = 7 log (u) 

\ {i- + J-a-l.s)^ + Y a+l,s) J 

15 In (2,2) signature the reality condition is (Y a ^ s (())* = F 0)S (l/£*). 

16 In this case the conditions (46) get modified to logli jS ~ — ^e~ — C s —D s for 9 — > —00 and logYi )S ~ 
— ^e — C s + D s for 9 —> +00. And similarly for the other Y functions. 
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For w 7^ 0, the 3x3 matrix A aa /(u) is invertible and we can multiply this relation by A^ 1 (w) 
to extract J-(l a , s )- For w = the matrix is not invertible and therefore when doing this 
operation we should allow for the constant zero modes in the final result. Finally, we rewrite 
the corresponding expression for F(l a ,s) i n position space. We obtain in this way the final 
set of integral equations 

log Y 2s = — m s y/2 cosh 9 — K 2 * a s — Ki* f3 s 

1 1 

logYi iS = -m s cosh9-C s --K 2 */3 s -K 1 *a s --K 3 *j s (47) 

log Y 3s = -m s cosh 9 + C s - \k 2 * (3 S - K x * a s + \k 3 * j s 



2 2 



where a,/3,7 are short-hand for 



a s = log jz ■ W1 , „ r , 7. = log 



;i + y 2 , b -i) (i + y 2 ,.+i) (i + Y l!S+1 ) (i + y 3 ,.-i 



:i + y 2 , s 



.2 



^ log (i + y, s -i) (i + iw) (i + y 3 , s -i) (i + y 3 ,,+i) ' (48) 

and the kernels read 

_ 1 v^ooshfl _ z tanh2 ^ 

1 27r cosh 6 1 ' 7T cosh 20 ' 3 7r 

The unusual appearance of a kernel which does not decay at infinity (K%) is a direct conse- 
quence of the singular behavior of A(w) at w — 0. 

Comparing the large asymptotics following from these equations with those predicted 
from the WKB analysis we see that the zero modes C s correspond precisely to the constants 
C s in (46) while the D s in the WKB asymptotics are given by D s = - J d9 r y s {9). 

A more straightforward exercise, compared with deriving the integral equations, is to 
check that they indeed yield the functional relations. To do so we simply compute the left 
hand side of the functional equations using the integral equations. When doing this we 
should use 

f ± = f(9±m/A T tO) 

in order not to touch the lines where the kernels K 2 and K 3 become singular. Then, simple 
identities such as + K± — K 2 = and + K 2 — 2Ki = 5(9) eliminate all the kernels 
in the right hand side of the integral equations and the functional equations are indeed 
reproduced. 

Up to now we have discussed the case where all masses are real and positive. To consider 
the case of complex masses m s = \m s \e liPs , we proceed in exactly the same way as described 
in section 3.5 for the AdS 3 TBA. That is, for small phases (p s the integral equations take the 
same form as in (47) with 

m s -> \m s \ , Y a!S (9) Y a , s (9 + i<p a ) , (9 - 9') (9 - 9' + i<p s - i Vs ,) , 
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where K stands for the three different kernels. At \ip s — (p s+ i\ = 7r/4, ir/2, 37r/4, ... we pick 
the poles from the appropriate kernels (see section 4.7 and appendix B for illustration). All 
in all, the Y's and therefore the area are continuous whereas the apparent jumps in the 
integral equations are just an issue of the choice of contour. 

4.5 Simple combinations of Y functions and s n+ i — > s\ monodromies 

When we normalize the solutions as in (35) it can happen that s n+ i is not equal to Si. Of 
course, they have to be proportional to each other. The proportionality constant is called 
a "formal monodrom". For n odd, this constant can be removed, by rescaling the solutions 
appropriately. For n even, there is some non-trivial gauge invariant information in this 
constant. In fact, this non-trivial information is a particular combination of Y functions 
which is particularly simple. 

This is most interesting for n = 4k so let us consider that case first. We introduce three 
combinations of Y functions 

B\ = (Y 1)n _ 5 Y 2)ri _ 6 Y lin _ 7 ) (Y 3n ^gY 2tn _ w Y 3n _n) 1 ( Yi i?1 _i3 Y 2]n _i4 Yl^-lb) ■ ■ ■ 

B 2 = (Y2 in _5Y" ln _ 6 y3 in _ 6 y 2 ,n-7)(^2,n-9^ / l,n-10^ / 3,n-10^2,n-ll) 1 • • • 

B3 = (Y^n-sl^.n-e^n-r) (Yi^-gl^/n-loYl^-ll) 1 (I3 ]n _i3l2,n-14 ^n-is) . . . 

They are the analogue of B((), equation (23), in the Ad S3 case. The Y-system alone, without 
recurring to the definition of the Y-functions, implies the following equations for B a : 



They have the form of a discrete Laplace equation for log(I? a ) with a non-diagonal metric. 
Given the expected analytic properties of this functions the solution to these equations is 



The value of wq can be found by computing the large 9 asymptotics in the definiton of the 
B's. We find that 



is the cycle around infinity and wq its complex conjugate. In our normalization these quan- 
tities turn out to be equal to B a = T atU _± hence we have not only derived the form (42) but 
we also computed w from the Y-system equations. 

For n = 4k + 2 we only have one simple combination which is 



B 2 B 2 _ B l B 3 _ B 3 B 
BiB A B 2 B 2 



1 . 




w = (m n _ 5 + \/2m n - 6 + m n _ 7 



) - (m„_ 9 + V2>, 



:m n _i + m n _n) + . . . 



n — 4 



k=l 



n 




(49) 
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This relation should be thought of as a gauge invariant definition of fi. The fact that fi is 
constant also follows directly from the integral equations. 

These relations imply some identities on the constants appearing in the WKB asymp- 
totics. By considering the small and large ( limit of (49) we find the expression for /i in 
terms of the C s 

n — 4 
2 

log// 2 = y^cwi 

k=l 

In addition, we find a relation on the D s of the form 

n — 4 
2 

= ^£ 2fe -i (50) 

k=l 

Note that we are viewing the C s as arbitrary constants, while the D s are determined by 
solving the integral equation. The results in this subsection can also be obtained by a direct 
analysis of the solutions at infinity, as is explained in appendix C. 

4.6 Area and free energy 

In this subsection we show that the area is related to the free energy of the TBA system. 
The final result is (58). The reader that is not interested in this derivation can jump to the 
final result. 

Note that the area is independent of (. However, our objective is to relate the cross 
ratios to the area. For this purpose, it is convenient to consider the small ( expansion of 
the Y functions. We have already encountered the small £ expansion when we looked at 
the asymptotic conditions. This expansion can be viewed as a WKB expansion where ( is 
playing the role of h. We take complex masses but with small enough phases so that the 
WKB approximations that we did before continue to be valid, with the same choice of cyl For 
our purpose, we need to expand these functions to first order in (. For small ( we diagonalize 
$ 2 = diag(xj), where Xi are roots of <iet($ 2 — x) = 0. In our case, we have x 4 — P(z) = 0. The 
Xi are different sheets of a Riemann surface over the z plane. Hitchin's equations imply that 
A s is also diagonal in the same basis. By a diagonal gauge transformation we can set A% — > 0. 
In general, A z is not diagonal in this basis. To order (°, only the diagonal components are 
relevant. Let us call these diagonal components by a\. Again, we can think of a as a one 
form on the Riemann surface. This is a closed form da = due to Hitchin's equations. 

When we evaluate the cross ratios we will be evaluating integrals of the form 

(i\U\i) = (i\Pe- fSi tr +A * dz+c ** d *\i) (51) 

Here U is the "evolution operator" taking us between two points in the z plane. The states 
|z) indicate that we follow a given branch, since changing branches is suppressed by an 
exponential amount. In other words, in the WKB approximation, we are following the 
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"ground state" and the excited states are integrated out and change the evolution (at order 
( and higher) of the ground state. Expanding to order ( we find that (51) has the form 



\og(i\U\i) 



c 



+ a l z dz + ((u l z dz + u\dz) + o(C 



(52) 



where u % is a certain one form. This one form is closed because the connection is flat. The u s 
components of this one form are simple, they only come from the diagonal components of $ s , 
u\ = $|\ The u z component comes from a term that involves the off diagonal components 
of A z , but we will not need the explicit expression here. 
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Once we form cross ratios, we find integrals over closed cycles on the Riemman surface 
of certain closed one form forms. Thus we obtain 



logK 



+ (h a + ( f u + o(C 



(53) 



where a and u are viewed now as forms on the Riemann surface (in the i th sheet they are 
a 1 or u l respectively). Since they were constructed from a flat connection, we see that the 
quantities in (53) do not depend on the precise choice of contour within its cohomology class. 
Thus, we can view them as conserved quantities of the Hitchin equation. In fact, we could 
expand to higher powers in ( and find extra conserved quantities, though (53) is enough for 
our purposes. 

In the gauge where we diagonalize we can rewrite the expression for the 



A = d 2 zTr[§ z <$> z ] = / d 2 zJ2 x i® 



xl dz A u\dz 



A 



xdz A u 



i 

-w 7tJ , 



xdz 



u 



(54) 



In the second line we expressed the integral as an integral over the Riemann surface. We then 
used a complete basis of cycles indexed by 7, and we used a formula that is valid when we 
integrate products of closed forms. Here w ia i is the inverse of the intersection form for the 
cycles (7,7/. This step is valid when the number of cycles is even (odd number of gluons), 
otherwise we will not find an inverse. In our particular problem, we will choose the basis of 
cycles to consist of all the WKB cycles associated to the Y a s functions. It is convenient to 
undo the shifts that define the Y functions and define Y atS functions which consist of various 
cross ratios evaluated at the same value of (, see (D.2) in appendix D. We then do the WKB 
expansion for a ( with phase e J7r//8 , so that all the Y a ^ s are associated to a corresponding cycle. 

It is useful to compute the small ( expansion of the Y a ^ s functions from the integral 
equations. We obtain 



log Y a>s 



T 



+ a a , s + C 



/dff 
e - e '\og(l + Y a , s (e>)) 
2ir 



(55) 
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If you are curious, the expression is u\ = ^2 k 



This is a familiar formula from second order 



■ik^i ( Xi ~x k ) ' 

perturbation theory (the Xi are the energy levels). From Hitchin's equations one can directly check that 
du = 0. 
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where M a ' s ' a > s , or M™ , is a certain matrix that results from the small ( expansion of the 
integral equation for Y a , s - By comparing (55) with (53) we can read off the values of § u. 
We obtain 

- j) Xdz = Z atS 

-i u = Z a>s + M a < s -> a '> s ' f f^e-* log(l + ? a>s (e')) 

When we insert these expressions in (54) we get two terms A = A periods + Af ree which are 
equal to 

Ap er iods g Wj^'ZyZ^i 

A free = -'-Z^w^M^'^' J ^-e- e '\og(l + Y a , s (e')) 

— Z a , s e- e log[l+Y a , s (e)} (56) 

a,s 

Here we have used that Z 7 w 7i7 /M 7 '' 7 " = — 4«Z 7 » for our case, due to the relation between 
the Z 1 (or m ajS ) with various values of a implicit in (46). The matrix of cycle intersections 
is given in appendix E, figure 17. We can take the average of (56) with a similar expression 
which we obtain if we did the large ( expansion to obtain 

/d6 
^[Z a , s e- e + Z ayS e e ] log[l + Y a , s (e)} (57) 

a,s 

A free = £ / f- K I cosh Q log [(1 + Y ljS ) (l + Y 3 , s ) (1 + F 2 , s )^] (9 + ia s ) (58) 

s J 

where e ias is the phase of m s . In the last equation we have written the final answer in terms 
of the Y functions (as opposed to the Y functions). We have also used the relation between 
m s and Z atS , implicit in (46) and explicit in appendix G.l, eqn. (94). The explicit value of 
Aperiods for in terms of the masses m s is given in appendix G.l. 

One can give an alternative derivation of this formula in the spirit of the one given in 
[15] for the n = 6 case. This alternative derivation starts with the observation that the 
area A can be viewed as the generating function of transformations that change (. This 
is the generating function when we define a Poisson bracket for the Hitchin system that 
makes $ 2 and $ 2 conjugate variables and similarly for A z and A s . Then one uses that the Y 
are variables whose Poisson brackets are computable and related to the cycle intersections. 
Finally one uses the integral equations as above to expand the Y functions for small (, to 
find the Poisson brackets of the quantities involved in this expansion. In this way one can 
check that the final expression for the area does indeed generate the desired transformation. 

Note that if we view the Hitchin system as arising from an M = 2 supersymmetric theory, 
then the Z 4 projection that we had would break M — 2 to J\f — 1 supersymmetry. This 
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M = 1 theory has a global symmetry for rotations of $ 2 and $ g in opposite directions. The 
area is then the D term potential (or momentum map) for this symmetry. This connection 
between a four dimensional supersymmetric theory and two dimensional quantum integrable 
models is in the spirit of [28]. It would be interesting to find the precise relation. Note 
however, that in our AdS 5 problem we do not have supersymmetry, we have integrability. 

4.7 The geometrical meaning of the Y-functions 

In the previous sections we saw how to determine the area of the minimal surfaces for a 
given choice of masses and chemical potentials in the F-system equations. To identify which 
polygons correspond to a given choice of masses and chemical potentials we must compute 
the space-time cross ratios (10) for these solutions. This can be done by evaluating the 
^-functions at special values of the spectral parameter as we explain in this section. In 
fact, the Y functions themselves are cross ratios, but written in terms of somewhat exotic 
variables, introduced in [29], as we explain in the next subsection, 4.7.1. In subsection 4.7.2 
we explain how to obtain the more conventional cross ratios defined in terms of distances 
between cusps. 

4.7.1 Relation between the Y functions and twistor cross ratios 

As we explained in the introduction, we can recover the form of the coset representative, 
g~ x by picking a set of independent solutions of the linear problem, see 4. In our case, 
this is a set of four solutions ip a of the linear problem. Each of these solutions is a spinor 
and a = 1,2,3,4 labels the solution number. We can orthonormalize them so that g^ 1 is 
a proper group element. These solutions then determine the shape of the string worldsheet 
in spacetime. Rotations of the a indices by a group element go corresponds to the SO(2,4) 
AdS isometries, g — > gog. A given solution, say i/; a , can be expanded in terms of four other 
solutions as 18 

1p a = X\s i+2 + r] l a S i+1 + ffcSi-! + fj* a 8i (59) 

where s i+2 is the big solution in the Stokes sector i, and Si±i are the intermediate solutions 
in Stokes sector %. As we a go to infinity within this Stokes sector s i+2 —> oo. Thus, this 
is the dominant solution and it determines the behavior of g~ l and eventually, also for the 
spacetime solution. We can calculate A„ as 

X\ = (Si-i, Si, Sj+i, 1p a ) = (Si, 1p a ) (60) 

Here X 1 is a spinor of the spacetime conformal group, which acts on the index a. 

Now, imagine that we want to evaluate inner products of the X 1 . In these inner products 
we contract the a indices of the four spinors with an e symbol. Notice that these are indices 
transforming under global conformal transformations. We obtain 

<A\ X j , X k , X 1 ) = e a > b > c '%XiX k X = (s h sj, s k , «,) 

18 It turns out that fj l = — A 4+1 . This follows by continuing this formula to the i + 1 Stokes sector and 
reexpressing s^_i in terms of the Sj appearing in (59) for that Stokes sector. 
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Here we have used that, by performing a local gauge transformation and a global transfor- 
mation, we can set the ip aA = <5 a ,a at some point on the worldsheet. 

Using this formula plus the identities (81), (82) we can show that 

( Si ,s i+1 , Sj ,s j+1 ) = (A\A l+1 ,A^ +1 > 
(s i ,s i ) = (s i _ 1 ,s i ,s i+1 ,s i ) = (X i ,X j ~ 1 ,X j ,X j+1 ) = (X\X j ) 

where X 1 = A* -1 A A 1 A A* +1 . Using these relations we can express the Y functions in terms 
of the A*, which are spacetime quantities. These "momentum twistors" can be introduced 
just from the knowledge of the position of the cusps [29]. In order to introduce them, we 
only need to know that we have a null sided polygon. When we introduce the A* from this 
latter point of view, the X 1 are defined up to an overall rescaling. In (60) we have picked 
a particular normalization. However, in the final expressions for the Y functions in terms 
of A, the overall normalization of each A* drops out, for the same reason that the overall 
normalization of the Si drops out. Each A* is associated to a null side of the polygon, and 
a pair of consecutive As determine the position of the cusp X % ab = A| a A^ , where X % ab are 
six coordinates defined up to a rescaling obeying X 2 = 0. Thus, they define a point on the 
boundary of AdS space. 

Note that the Y"i iS and Y 3tS only differ by s -h- s, or A -H- A, see appendix D. This operation 
is simply target space parity. 



4.7.2 Traditional cross ratios from the Y functions 



In this subsection, we explain how to obtain traditional cross ratios from the Y function. By 
a "traditional" cross ratio we mean one constructed from physical distances as in (10). These 
can be introduced via 



1 2,s 1 2 y 



6»=j7rr/4 



where we combined the definition of V-function in terms of the T-functions with Hirota 
equation (39). This ratio has been constructed so that it only involves the functions T2 jS . 
These are determinants of four small solutions of the form (sjSj + x%'%+i)- We recall that the 
physical cross ratios are ratios of four such quantities (10). For example, 

jj[Q\ _ ( s -k, s -fc+l) Sk, Sfc+l) (S-fc-l, Sfc_i, Sfc) _ x^ fc fc x 2 k _ l k _ l 

2k— 2 7 W \" 2 2 ' V / 

\ s -k-l, S-fc, Sfc, Sk+l) \ s -k, S-k+li s k-l, Sk) X_ fc _ l fc X_ fej j fc _ 1 

where Xjj = Xj — Xj, see figure 8. If we consider U^_ 2 then we will shift the position of the 
small solutions by p sectors, i.e. we will get a cross ratio involving the cusps {x_ fc _ 1+p , x_ fc+p } 
and {xfc_i +p , Xfc+p} of the polygon. We see that the index 2k — 2 is the number of cusps 
between x_fc +p and Xfc +P _i (counted from the side of the cusp x ). Similarly, the cross ratio 
involving sides {x_ fc _ 1+p , x_ fe+p } and {x fe+p , x fc+1+p }, which are separated by an odd number 

of cusps, are given by uf k p ^\ 
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X-k-l,k-l 




Figure 8: The cross ratio U* 




-fc,fe x -fe-i,fc-i 



A generic cross ratio involving four non-consecutive cusps of the polygon can easily be 
constructed by multiplying appropriate factors of U\ m \ e.g. 



We see that once we solved the integral equations it suffices to evaluate the ^-functions Y 2 , s 
at some specific values of 9 to read off the corresponding cross ratios. 

A subtlety of pragmatic nature is the following. When solving the integral equation we 
usually do it numerically by iterating the integral equations. At the end of the iteration 
cycle we are left with the functions Y ajS (9) in the real axis which is precisely what we need 
to compute the free energy (58). On the other hand, to get the physical cross ratios we will 
generically need to evaluate Y 2>s at some imaginary values. Using the integral equations, 
these cross ratios can be written in terms of (integrals of) the F-functions evaluated in the 
real axis. When continuing the integral equations (47) out of the physical strip we will 
have to pick some extra pole contributions from the several kernels as we cross the lines 
Im(6>) = ±27r/4, ±?7r/2 etc. 

To make this procedure clear let us illustrate how to compute Y 2 ^(in /2). First we notice 
that Y a ^ s (in/ 4) can still be expressed in terms of integrals of the Y- functions over the real 
axis using the original undeformed equations (47). The only thing we must be careful about 
is that the kernels K 2 and K3 have a pole singularity for 9 = in/ 4 so that we should either 
interpret the integration contour to go slightly above the real axis or equivalently we should 
plug 9 = z7r/4 — iO in the right hand side of these equations. To reach values of 9 with an 
even larger imaginary part (such at 9 = m /2) we can simply pick the pole singularities of 
these to kernels. For example 

log (Y 2jS ) (9) = -m s V2cosh(9) - K 2 * a s - K x -k /3 S + a s {9 - in/ 4) , tt/4 < lm(0) < tt/2 

We can now evaluate the right hand side of this equation at in/2 — iO to compute Y 2yS (iir/2). 
The last term becomes a s (iir/4:) and contains a bunch of ^-functions evaluated at 9 = in /A. 
We already explained how to get those in terms of integrals of the ^-functions in the real 
axis using the original equations. Hence we are done. 




IK 



1=1 
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There is an even more efficient way of computing these cross ratios which goes as follows. 
Suppose we need some Y-function Y atS (6) where 9 is outside the physical strip, i.e. |lm(0)| > 
7r/4. For concreteness let us suppose Im(#) < — ir/4. Then, by repeatedly using the functional 
equations (45) as 

Y _ (1 ^a,m+l)(l + ^4-q,ro-l) 

we can express Y ajS (9) purely in terms of Y-functions inside the physical strip. For those we 
can use the original integral equations (47) as explained above. 



4.7.3 ( symmetry 

We explained how to get the physical cross ratios by evaluating the Y-functions at 9 = 
0, ±i7r/4, ±i7r/2 etc. Actually we can construct a family of polygons parametrized by a 
complex number 9q which will all have the same area! These polygons are obtained by 
evaluating the Y-functions at 9 = 9q,9q ± 27r/4, #o ± m /2 etc. The reason why we can 
evaluate the y-functions around any # is that flat sections of the linear problem can be 
assembled into a physical solution for general ( = e e ° at not only for ( = l. 19 For purely 
imaginary 9q the new polygons are real, for generic 9q they are complexified. Thus we can 
view changes of ( as a symmetry of the problem. Namely, we have a one parameter family of 
polygons labeled by ( which have the same area. It would be very interesting to understand 
this symmetry in greater detail and specially to see if/how it manifests itself at the quantum 
level. 

Note that this is intimately related to the fact that changes in ( are generated by the area 
when one defines a certain Poisson bracket that is natural in the theory of Hitchin systems 
[15]. This Poisson bracket is different from the one one gets from the sigma model [35]. 



4.8 High temperature limit 

In this section we will focus on a particular kinematical regime. From the Y-system point 
of view we want to consider the limit when the ^-functions are approximately constant (in 
some large region of 9). To find the constant values of the Y-functions we solve the Y-system 
equations (45) and (49) dropping the ( dependence, Y^ s — > Y ajS , 

Ya,mY&—a;m (1 Yi,m+l)(l Yt— l) (62) 

Yi+l,mYi— l,m (1 Yt+l,m)(l Yi— l,m) 

For this approximation to be valid, the sources in the integral equations must become in- 
dependent of the spectral parameter which means m s — > 0. In the TBA context this arises 
in the high temperature limit, so will adopt that terminology here. Of course, in the AdS 
problem there is no temperature since we are just solving classical equations. The m s — > 
condition is however not sufficient. Namely we do not have the freedom to chose the values of 

19 Because £ ^ 1 can be absorbed into a redefinition of P(z) — > £ 4 P(z) together with a ( dependent global 
gauge transformation. 
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the constants C s if we want the solution to the Y-system to be given by constant Y-functions. 
Instead these constants are found from the Y-system. 20 To derive this, note that the ratio 
of the equations (62) for a = 1 and a = 3 implies that 7 S = where 7 S is defined in (48). 
Thus, from equations (47), we see that 



2C S _ YM. 

Yi,s • 

We will first consider the richer case of an even number of gluons. The n odd case is 
discussed afterwards. The equations (62) admit a one parameter family of solutions. We can 
parameterize this family in terms of 

n — 4 

For a fixed value of fi, (62) admits a discreet set of solutions. At n — 1 there is a unique 
solution such that Y ajS > for any a = 1,2,3 and s = 1, ... ,n — 5. We will consider a 
solution valid for arbitrary [i and continuously connected to that unique solution. This one 
parameter family describes a family of regular polygons with a Z n symmetry. This is shown 
in more detail in appendix I. 

The polygon can be described in terms of the n twistors 

A fc = {(1 + a k ) (1 + b k )e^, (1 + a k )e^ , (1 + b k )e~^, e~~ } (64) 
where fi = e 1 ^ and 

. 27T (f) . , 71 (j) 

a k = (— 1) tan — tan— , b k = (— 1) tan — tan — . 

n n n n 

Then the constant solution to the Y-system is simply obtained by plugging the twistor \ k 
in place of s k in the relations (84) in appendix D. 

When = we have a k = b k = and the solution simplifies to 

sin(^(4-a + s))sin(^(a + s)) , 
sin g(4 -a)) sin (0=) 

Actually, in this case we have Y a>s = Y 4 _ a s and hence the Y-system reduces to the standard 
Y-system and the solution (65) can be found in the literature [30]. In fact, this is a regular 
polygon that can be embedded in A0IS4. Another limit where we find significant simplification 
is the limit where <fi — » ix(n — 4)/2. In this limit we find Yi^fc+i = Ys^+i = —1, Yl,2& = 
Ys,2fc = 0, Y 2i2 fc+i -> 00 and 

k + 2\ . f k\ , 2 /7r\ A n 



Y 2i2fc sin ( 7T— — 1 sin ( vr^ J /sin (^-J , fi = - . (66) 



20 In this sense the nomenclature high temperature is slightly abusive. The genuine high temperature limit 
would correspond to m s — > 0, with the chemical potentials C s arbitrary. 
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This corresponds to a regular polygon that can be embedded in AdS 3 . In fact, the curious 
pattern of zeros and infinities for the F-functions is a generic feature of the AdS$ limit. In 
the next section 4.9, this is precisely how the F-system in the AdS$ strip reduces to the 
^-system in the AdS^ line. As we move between = and = ! ^7r we interpolate 
between these two cases. 

It turns out that the free energy can be computed exactly in this limit, as shown for 
instance in [31], and reads 

F - = E H(e (2 ~ a)Cs Ya,s) log(l + Y aiS ) + 2Li 2 (-Y a , s )] . (67) 

a,s 

Plugging the analytic expressions above into this expression we find the remarkably simple 
result 

F„ = £((n-4)(n-5)-^) . (68) 

As already mentioned, for = we recover the regular polygons of second class that can be 
embedded into AdS±. One can actually see that the free energy for = exactly reproduces 
the numerical prediction (98). Similarly, for = ! ^vr we can see that (68) exactly reproduces 
the AdS 3 result (99). 

In addition to the solutions described above, there can be discrete families of solutions. 
For instance, let us focus in the n-odd case. In appendix I we have constructed a set of 
solutions parametrized by the number of sides n and an extra integer r, with r = 2, (n — 
l)/2. The spinors characterizing such solutions are 

^ = ( e w(r+l)^ e t ff (r-l)^ e -tir(r-l)^ e -tir(r+l)*\ 

where we have reorganized the expression appearing in the appendix. From these spinors, 
one can easily compute the cross-ratios Y as . As these polygons can be embedded into AdS±, 
we expect /i 2 = 1 and indeed we find ^ = 1. In addition, one can explicitly check that this 
cross-ratios give a solution of the Y — system equations. This constitutes a non trivial check 
of the Y — system for the case in which n is odd. 

Let us comment on an interesting feature of (68). Note that A in (64) is a periodic 
function of with period mr. This means that the cross-ratios describing the polygon are 
periodic functions of 0. On the other hand we see that (68) is not periodic! We thus have a 
family of solution all ending on the same polygon. If the right prescription is to sum over all 
these surfaces, then the full result will have no monodromy. In that sum, one solution will 
dominate while the others are non-perturbative corrections (in l/y/\). On the other hand, 
in terms of amplitudes we expect non-trivial analytic continuation properties. For example, 
it is well known that amplitudes have interesting monodromies as we analytically continue 
the external parameters. This is an essential tool in weak coupling computations, see [32, 33] 
for example. If the right prescription is then not to sum over all these surfaces, then it would 
be interesting to study this particular kinematic configuration and understand in a better 
way the physical meaning of the monodromy in (68). 

Such non-trivial monodromies for the free energy are a common occurrence in TBA 
systems. After an analytic continuation in the parameters we do not end up with the 
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ground state, but we end up in an excited state [34]. In fact, the high temperature TBA 
typically corresponds to a CFT. Then the chemical potential simply translates into a winding 
condition for the scalar that bosonizes the corresponding U(l) current and thus giving the (p 2 
contribution to the energy, see e.g. [36, 37]. This, in particular, shows that in our problem 
some excited states will typically appear as we analytically continue the parameter [36]. 
These are related to poles in the physical region and need to be taken into account according 
to a well understood procedure [34]. Namely, we need to deform the contours and pick up 
the corresponding poles, etc. 

4.8.1 High temperature limit of AdS 3 Y— system 

We can consider the high temperature limit of the TBA equations corresponding to scattering 
amplitudes on AdS% for n = 2n gluons. In this case the F-system equations (18) becomes 
simply Y 2 = (1 + Y s+ i)(l + and the constant solution reads 




These are precisely the F-functions found in the previous section in equation (66). This 
is not surprising since we had already anticipated that for <p — > (n — 4)ir/2 the polygons 
described by the previous solution become AdSs solutions. The AdS^ free energy can be 
then evaluated, exactly as in (67), and we obtain 

f« = 4E wn) + F *) + 2L ^(-Ys)) = ^-g)^- 4 ) (69) 

Z7r * — ' LZn 

s 

where we have reinstated n = 2n. This result is however not the same as (68) for = (n — 
4)7r/2 even though we are describing the same solution. This is actually not a contradiction: 
as explained in appendix H the regularization of the area in the high temperature limit 
amounts to subtracting the area of n — 4 regular pentagons in the AdS^ case and the area of 
n/2 — 2 regular hexagons in the AdS$ case. Taking this into account, the difference between 
two free energies is precisely as expected! 

4.9 AdSi and AdS^ reductions 

Minimal surfaces that can be embedded in an AdS^ or AdS% subspaces of AdS§ are more 
restricted and as a result, the problem simplifies. The reduction of the AdS§ flat connection 
was done in [15]. In this section, we will consider the implications of this simplification to 
the Y-system. 

The worldsheet theory describing strings moving in an AdS^ subspace is obtained from 
the parent AdS^ by an additional projection. This projection relates A to A 1 via a gauge 
transformation F . That gauge transformation therefore relates solutions to the problem 
with solutions to the inverse problem. More precisely we have Sj(£) = Fsi(Q, where in 
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our normalization det(F) = 1. As a result, Yi jS = Y 3 ^ s (and hence /i 2 = 1). The Y-system 
equations can then be written in the form 



^2,m^2,m (1 + ^2,m+l) (1 + ^2,m-l) ^\,rrX\,m (1 + Yl,m+l) (1 + Y\ } 



m—1 j 



A solution in AfiS^ that can be embedded in AdS% must have even number of gluons. 
The linear problem splits into two decoupled problem denoted by left and right problems in 
[17]. In an appropriate gauge we can write 

( S R,k \ ( 

S2k — \ n , S2fc+1 — 

where sl and sr are the small solution of the left and right Ad S3 problems respectively. 
Because of this factorization the T-functions can be dramatically simplified. We choose 
a normalization where the AdS 3 solutions obey (s^i, = 1, and the same for the 

right. 21 The left and the right problems are related by a rotation in the spectral parameter 
(sR, a , SR : b) = (sL,a, We can use this relation to translate all inner products into inner 

products of the left problem. We also recall the definition, T k = (sl,q, SL,k+i)^~ 2k ~ 2 \ of the 
Hirota functions for the AdS% problem in section 3.2. We then find that the the Hirota 
variables T a>s of the AdS§ problem have the form 

Tl,2fc+1 (C) = , Tq jS = T^ s = —1 , 



T 1)2k (0 = -(s Lfi ,s LM1 )^ = -T, 



k ' 



r 2i2fc (C) = -(^ 1 ^ 1 > [ -*- 11 (« IflI ^i) [ -*- 11 = -rW 1 
T 2 , 2fc+1 (C) = (^,o,^> [ - 2fc - 2 k^ 

and of course T 3tS (() = Ti iS (£). 

The AdS§ Hirota equations (39) becomes identically satisfied in all nodes except s = 2k 
and a = 2. For these, it becomes 22 

{T k f (T [2] T [ k 2] - T k ^T k+1 - l) = , 

Inside the parentheses we recognize Hirota equation (16) in AdS%. Recall that in the AdS 3 
treatment supercripts denoted shifts in the spectral parameter which were twice as large 
compared to the ones we are using now, e.g. 



{f^ ')/iere {f ) 



AdSs section 



The 3n — 15 nodes in the AdS$ strip reduces to the n/2 — 3 nodes in the AdS% line in a very 
funny way which we represent in figure 9. Basically the only non-trivial functions become 

r 2| 

^22fc = Tk+iTk-i = Y k and these obey the AdS 3 Y-system equation (18). 



21 Note that with this choice, the AdS§ solutions obey (sj, Sj+i, Si+2, Si+3) 
22 We actually obtain this expression at e I7r / 2 £. 
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1 2 n/2-3 




Figure 9: In the AdS^ reduction most of the 3n — 15 functions Y a ^ s become trivial (either zero, 
one or infinity as indicated in the corresponding node in the figure). The n/2 — 3 nontrivial 
functions Y^fc = Yk are indicated by red circles. The AdS$ F-system in the original strip 
(with shifts of m/A) imply the F-system equation for the line of red nodes (with spectral 
parameter shifts of mt/2). 

5 Conclusions 

In this paper we have given a way to compute the area of minimal surfaces that end on a 
null polygon at the boundary of AdS$. The method uses the integrability 23 of the classical 
equations in an essential fashion. 

We have used the map between the integrable sigma model and a Hitchin system. The 
Hitchin system is an SU (4), or SU (2, 2), or SL(A) system depending on the signature. More 
precisely, it is a certain Z 4 projection of this system. Alternatively, we can simply say that 
we have chosen a specific form for the one parameter family of flat connections. This family 
is parametrized by a spectral parameter ( = e e . 

For this problem the worldsheet is the complex plane and we have an irregular singularity 
at infinity. This means that there are Stokes sectors as we approach infinity. Each of 
these Stokes sectors has an associated small solution Sj and a large solution. The large 
solution determines a four-spinor Aj, which specifies the direction in which the large solution 
is pointing. These spinors are associated to the sides of the polygon and are the same as the 
momentum twistors introduced by Hodges [29]. Alternatively, we can say that consecutive 
Stokes sectors determine a cusp or a vertex of the polygon. Using this cusp positions, or 
the momentum twistors A, we can construct cross ratios. We can introduce a family of cross 
ratios depending on the spectral parameter (. If £ = 1 we recover the physical cross ratios 
of the original problem. 

A particular set of cross ratios, denoted by the Y functions Y a>s , obeys a set of functional 
equations, or Y system equations (45). The number of Y functions is 3(n — 5), which 
is the same as the number of independent cross ratios of the problem. These functional 
equations, together with some input regarding their asymptotic properties, imply a set of 

23 The Yangian charges - responsible for integrability - are also visible at weak coupling [39]. Still, to this 
moment it is not clear how exploit integrability efficiently in that regime. At strong coupling the connection 
between conformal and dual conformal symmetry and integrability was worked out in [40]. 
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integral equations for the Y a ^ s functions (47). These integral equations involve 3(n — 5) real 
parameters. These parameters were not present in the Y system function equations but 
they appeared in the specification of the boundary conditions for the Y functions. The 
solution to these equations relates these parameters to the physical cross ratios. Roughly 
speaking these parameters are related to the values of Y at 9 = ±oo, while the physical cross 
ratios are Y{9 = 1). These functional equations have the form of Thermodynamic Bethe 
Ansatz equations for a certain quantum system which is not associated in any obvious way 
to our initial problem, which was a classical problem. Similar relations were observed in [38]. 
Moreover, the area is given by the free energy of the TBA system. In practice, the TBA 
equations can be solved numerically. As in [15], we have been able to to solve the equations 
in a specific "high temperature" regime. This gives a one parameter family of solutions (for 
n even) . These describe a family of regular polygonal contours at the AdS boundary. This is 
a one dimensional line in the 3(n — 5) dimensional space of cross ratios. The answer for the 
area is surprisingly simple, it is just A oc <p 2 , where <fi parametrizes the family and it appears 
as a chemical potential of the TBA system. The cross ratios are periodic functions of <fi. This 
simple form is expected from the TBA perspective, since it is describing the UV limit of the 
2d quantum integrable theory which is a CFT, and <j) appears as a chemical potential [37]. 
One interesting aspect of this family of solutions is that, since we know it analytically, we can 
analytically continue in the space of parameters and find interesting monodromies. Namely, 
the cross ratios come back to themselves while the area does not. It would be interesting to 
study them further. One thing we can say, is that this shows that we will get TBA equations 
for excited states as we do analytic continuations. By now, this is a familiar phenomenon 
[34]. Thus the full problem involves not only the TBA ground state but also some excited 
states. 

There are several interesting problems for the future. One of them is to take the large 
n limit in order to obtain arbitrary (spacelike or timelike) contours. This would effectively 
solve the strong coupling form of the loop equations [41]. 

Of course, the most interesting open problem is the extension of this to the full quantum 
theory. This will probably require, as a first step, the knowledge of the classical solutions 
for the full AdS& x S 5 sigma model. 

We have emphasized that we get the physical values of the cross ratios by evaluating 
the Y functions at ( = 1. However, we also get equally nice, but different, physical values 
by taking £ = e lLp . By varying (p we move in the space of cross ratios. All these values of 
the cross ratios have the same area!. Thus, changing ( corresponds to a symmetry of the 
problem. Namely, by changing ( we change the cross ratios in a way that does not change 
the area. Other values of £, with |£| ^ 1 correspond to generically complex values of the 
cross ratios and represents an analytic continuation of the problem, an analytic continuation 
that keeps the area fixed. Recall also that, with a certain definition of Poisson brackets, the 
area is precisely the generating function for this symmetry [15]. 

One curious observation is the following. We have the formula for the amplitude as: 
Amplitude = e~^^ Area \ Since the area is the free energy, this formula looks like we are 
computing the partition function of the system on a torus, where one of the sides has length 
proportional to vA- For large vA only the ground state contributes, which is what we 



40 



Figure 10: Left: Regularized area as a function of two real masses mi and m 2 for the decagon 
with n = 10. When the masses are very large the V-functions are exponentially suppressed 
and the free energy vanishes. When they are very small the free energy tends to the analytic 
prediction F = 7r/5 ~ 0.63 from the high temperature limit. Right: For those values of m\ 
and rri2 the cross ratios Yi(0) and 3^(0) cover a diamond shaped region. The regularized 
area is plotted as function of these cross ratios. 

computed. The overall sign is not quite right for this interpretation. It is nevertheless 
suggestive. 
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A Numerics 

In this section we explain how to implement equations (25) numerically in Mathematica in a 
very simple way (the code is quite similar to the one used in [42]). The algorithm is trivial, 
we simply iterate the integral equation plugging the V-functions at iteration k — 1 in the 
right hand side of (25) and reading from the left hand side their values at the fc-th iteration. 
We start by defining the kernel appearing in the integral equations, 
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K[x_]=l/(2 Pi Cosh[x]); 



and specify how many gluons we want to consider. We do that by introducing a list of 
masses which appear in the integral equations. For example, we set the following numerical 
values of the masses 

m={l . ,2.}; M=Length [m] ; 

would correspond to M = 2 nodes, i.e. to a polygon with 2M + 6 = 10 sides. We also 
introduce a cut off for the several integrals and set the number of iterations, 24 

cut = ArcCosh[8Log[10]/Min[m]] ; ni = 8; 

At each iteration step we compute the new values of the ^-functions at a discrete set of 
points and construct the new ^-function as the function which interpolates through these 
points. For that we need an interpolating function and a function to preform the several 
integrals 25 , 

F[S_] := FunctionInterpolation[S,{x,-cut,cut}, InterpolationPoints->100] 
int [S_] : = NIntegrate [S , {y , -cut , cut}] / ; NumberQ [S/ . y-> . 1] 

This is all the machinery we need. The iterations are simply implemented by 

Y[a_,k_] := (0#&)/;a==0| |a==M+l 
Y[i_,l] := Exp [-m[[i]] Cosh [#]]& 
Y[i_,k_] := Y[i,k] = 

F [Exp [-m [ [i] ] Cosh [x] +int [K [x-y] Log [ ( 1+Y [i-1 , k-1] [y] ) (1+Y [i+1 , k-1] [y] ) ] ] ] ] 

In the first line we set Y = Ym+i = 0, in the second line we set the seed values for the 
iteration procedure to be simply the asymptotic expressions for the V-functions and in the 
last command we compute the ^-functions at the iteration step k using the right hand side 
of the integral equations with the y-functions at step k — 1. We can now compute the free 
energy from (33) and see how it converges to a particular value as we iterate: 

energy=Table [int [Sum [Cosh [y] m [ [i] ] Log [1+Y [i , s] [y] ] /Pi/2 , {i ,M}] ] , {s , ni}] 
ListPlot [energy] 

For the masses chosen above we see that the free energy converges to 0.27. Finally, if we 
want to read the cross-ratios of the polygon whose area we have just computed, we simply 
need to evaluate the F-functions at 9 = 0, 

24 This cut-off is chosen so that for the smallest mass we have e - mcosh ( Cut ) — 10~ 8 . We could do something 
fancier (and more efficient of course) and introduce a cut-off which would be different for different integrals. 

25 The last command in the integration function ensures that we are really integrating something which is 
numerical once the integration variable y takes some random numerical value. We chose 0.1 but anything 
would work. 
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Table [Y[i,ni] [0] ,{i,M}] 



To compute the the free energy again for a different choice of masses we should first run 
Clear [Y] . 

The results of the numerics are shown for several values of mi and m2 in figure 10. Here 
we have considered real values of mi, m 2 . One can similarly do computations for complex 
values of mi, m 2 , using the kernels after (26), and the discussion in appendix B, if necessary. 



B Wall crossing and TBA 

In this appendix we describe the pattern under which the AdSs integral equation (26) changes 
as the phases of the masses (p s = arg(m s ) are deformed beyond the region where \(p s — (p s +i \ < 
7r/2. The AdSz integral equations (47) follows an analogous pattern. In the context of the 
works [26, 25] this corresponds to the wall crossing phenomenon. 

As explained in section 3.5, as long as \tp s — <£ s +i\ < 7r/2, the integral equation reads 

\ogY s (6) = -Kl cosh# + K SyS+l *log(l + Y s+ i) + K s ^ t *log(l + Y s - X ) (70) 
where Y s {6) = Y s (9 + iip s ) and 



K s . s , * log(l + Y s i 



r \o g (i + y s ,(6>)) d0l 

J 2tc cosh(0 — 6' + i(p ss ') 



and ip s ,s±i = fs — f s ±i- Let us consider the simplest possible situation where <^i,2 crosses 
the value +7r/2 (and hence tf2,i crosses — vr/2) while all other |<^ S)S+ i| < n/2. Any other 
possibility can be treated similarly. 

When (p x> 2 crosses the value +tt/2 we pick the pole contribution from the kernel and the 
TBA equation for s — 1 becomes 

hgYi(6) = -|m 1 |cosh^ + ir li2 *log(l + F 2 ) + log(l + y 2 , (^)) , (71) 

where 

Y^9) = Y 2 (9 + i<p lfi -m/2) 

is a "new" Y-function which arose in this process. Such Y-functions arise in pairs. Indeed, 
at the same time (f2,i crosses the the point — vr/2, so that the equation for Y% becomes 

log? 2 (#) = -|m 2 | coshtf + ^2,s ^log(l + Y s ) + log(l + Y{{9)) , (72) 

8={1,3} 

where 

Yl(6)=Y 1 (6 + i<p 2il + tn/2) . 
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These equations are correct for 7r/2 < </? lj2 < ir. We evaluate (71) and (72) at 9 =F 2(^1,2 — 
7t/2) ±z0 to get 

logY/^) = -|mi|cosh(0-i^ 1)2 + i7r/2) + ^i )2 *log(l+y 2 )+log(l + r 2 (0)) (73) 
logY 2 '(0) = -|m 2 |cosh(0 + i^-i7r/2)+ K 2 , s *log(l + Y.) + log(l + ?i(0)) , 

s={l,3} 

where X's are the kernels shifted accordingly. 26 

These equations, together with (71), (72) and (70) for s > 2 constitute a closed set of 
n — 1 equations for the functions Yi,Y 2 , . . . , Y/, Yj'. If more phases become large we get 
in a similar way more extra functions Y,'. 

As usual in integrable models, nothing singular happens when the phases cross the values 
of ±7r/2 and new singularities hit the integration contour. After all the mechanism of picking 
the poles is precisely designed to keep everything smooth!. In particular there is no reason 
to change the expression of the area (33) which therefore keeps its form. Notice in particular 
that the tilded Y-functions are designed to be always exponentially suppressed as 9 — > ±00. 

Wall crossing in the AdS$ system can be treated exactly in the same way. 



B.l Relation to Wall Crossing in [26] 

In this same circumstance, the authors of [26] obtained only one extra function rather than 
two functions as we obtained above. This single extra function is associated to a new cycle 
on the Riemann surface (in the WKB approximation) which corresponds 2 ' to 71+2 = 71 +72- 
Thus, Z 1+2 — Z\ + Z 2) or m 1+2 = m 1 + im 2 = |^i+2|e*' ,0l+2 , where the rrii are the complex 
m's. In the language of [26] a new hypermultiplet has appeared. 

In fact, we can define a new set of functions 



i + Y£ l + yj 

yn = Y i( e + m+2 - m)Y 2 (9^+ upi+2 - m - W 2 ) ^ 

1 + Yi(9 + i(pi +2 — i<Pi) +Y 2 (9 + i(fi+2 — — in/2) 

The first two variables are obviously designed to absorb the log(l + Y') terms in the right 
hand side of (71) (72). These relations look simpler when expressed in terms of the Y 
functions (without the tildes) 28 



+ y- ' 2 1 + y+ ' 1+2 1 + y 1 + y 2 



26 We can absorb the log(l + Y" s )'s in the right hand side of these equations by flipping the sign of iO in 
K12 and ^2,1 ■ This iO prescription is necessary to define these shifted kernels precisely. Otherwise there is 
singularity along the integration contour. 

27 Here we will use 1 + 2 as an index, hopefully, this will not cause confusion. 

28 Of course, they are even simpler in terms of the Y s functions defined in (28), since the shifts disappear 
[26]. 
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where the ± index is the usual shift by ±z7r/2. The function Y{\ 2 was defined so that 

(1 + F) = (l + r«)(l + F« +2 ) (75) 

(i + f 2 ) = (i + r 2 n )(i + Yi" + + 2 ) 

which transforms (71) (72) into equations that look more like the equations in [26] 
log Y? = -K| cosh + K 1>2 * log(l + F 2 n ) + +2 * log(l + y^ 2 ) 
log F™ = - \m x | cosh + K ^ * log(l + Y s n ) + K 2 ,i + 2 * log(l + Y£ 2 ) , 



s=l,3 



logF 3 = -\m 3 \ coshtf + ^ K 3 , s *log(l + Y?) + tf+ 1+2 *log(l + Y? +2 ) 

s=2,4 

where F™ = Y s for s > 2. The rest of the equations remains the same. Finally, the 
new equation for YJ\ 2 follows from considering the sum of (71) and (72) evaluated at the 
appropriate values, which are the ones in the numerator of (74). One of the kernels gives a 
delta function and produces a factor of log(l + Y™ +2 ). This combines with other terms to 
give 

log Y? +2 {6) = -|mi + im 2 \ coshtf + K w * log(l + F") + K w * log(l + Y 2 n ) 

+^r +2j 3*iog(i + y 3 ) (76) 

In the previous subsection, the equation for the free energy did not change, it was still 
given by Y\ and Y 2 , with no appearance of two extra functions Y(, Y 2 . Here, however, due 
to (75), we have a change in the expression of the free energy to 

/(1ft r ~ ~ ~ 

— cosh^ log(l + F") + |m 2 | log(l + Y 2 ) + |m a + im 2 \ log(l + F" 2 ) 
27T L 

plus the usual terms for s > 2. From our point of view, it is not clear whether this way of 
writing the integral equation has any advantage with respect to the one we wrote above. 



C Asymptotic behavior of the solutions at large z 

In this appendix we prove some formulas that are necessary to derive the AdS^ F-system 
equations. 

We will start by defining the small solutions and in a certain normalization which 
simplifies some of the formulas. Of course, the final F system involves cross ratios and 
is independent of such normalization details. Up to a normalization, Sj is a solution of 
(d + A(())s = that is small in Stokes sector i. Similarly, Sj is a solution of (d — *4*(C))s = 
which is small in Stokes sector i. Our first goal is to set a normalization of all these solutions 
so that (35) are valid. 
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It is convenient to introduce the w "plane" via dw = P 1 l 4 dz. For large z, or large w we 
cover the w plane n/A times as we go around the z plane, since for large z we have w oc z n ' 4 , 
where n is the number of gluons. 

Our boundary conditions for the connection are such that, for large w the connection is 
d + A(C) ~ d + diag(l, —i, —1, + diag(l, i, —1, —i)(dw (77) 

We can now choose a basis of approximate solutions of the form 

V>i = e" (M,/c+,I)C) (l, 0,0,0)* 

^ 2 = e -(-wc+™c)( 0; lj0) o)* 

^ 3 = e- { - w/<: -* o {0, 0,1,0)* (78) 
^ 4 = e -(i«/C-i*O( 0> 0,0,1)* 

For zero phase of (, the above solutions are the small solutions in consecutive sectors 
starting from the one centered on the positive real axis, followed by the one centered on the 
positive imaginary axis, and so on. 

We normalize the small solutions Si by saying that Sj = ip a with a = imod(4) in the 
corresponding Stokes sector. 29 Similarly, we choose a basis for the bar solutions 

= e (w/c+,i}C) (l, 0,0,0)* 
^ 2 = e ( -™ /c+ ™ c) (0, 1,0,0)* 

t/j 3 = e (- w K-™ c \o, 0,1,0)* 

^ 4 = e (WC-i<BC)( 0j 0,0,1)* 
The normalization of Sj is fixed by setting 

Si = tp a+2 , i = a + 2 mod(4) 
in the corresponding Stokes sector. 
We can now check that 

i^i, Sj+l, Sj + 2; s i+3) = ( — 1)* +1 

by evaluating this expression at infinity in sectors % + 1 or % + 2 where the asymptotics (78) 
of all these small solutions are still reliable. Similarly, 

Si = (-l) i+ Vi A Si As i+l (79) 



s k+1 (C) = (C) T s k (i() 

s k+1 (C) = (Cy's^Q ; C- 1 



( 1 \ 

10 

1 

\ 1 / 



29 In principle, we could analyze Si in the Stokes sector j, but then Si will be, in general, an arbitrary 
linear combination of (78). 
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This is not exactly as in the main text, so we perform a further redefinition of the small 
solutions 

$Ak+a { — ^) k SAk+a , SAk+a+2 ~^ ( — ^) k $Ak+a+2 , a = 1, 2, 3, 4 , (80) 

This transforms the relations (79) into the ones in the main text (35). It also transforms the 
matrix C into the matrix C in the main text, see (34), with det C = 1. 

Using epsilon symbol identities one can check (easily with Mathematica) that 

(s k , Sfc+i, sj, s j+ i) = (sk, Sfe+i, sj, s j+ i) (81) 
(Sfc> Sfe+l, s m ) = (Sfe+1, s m _x, s m , s m +i) (82) 
This identities, together with (35), imply that 

(s fc ,s fc+ i,s i ,s i+ i)(C) = (s fe _i,sjfc,s i _ 1 ,Sj) (e l7r/2 C) (83) 

(sj,Sjfc +1 )(C) = (sj,s fc ,s fc+ i,s fc+2 )(C) = (s i ,s j _i,s j _2,s fe )(e m/2 C) = (sj-i, s fc )(e i7r/2 C) 

These two relations are the important identities which we use in the main text. In deriving 
this relation one needs to use that if U is any matrix with unit determinant, then the product 
of four solutions {Usi} is the same as the product of the four solutions {sj}. 

Let us conclude with some remarks about the monodromy when % — > % + n. Notice that if 
n 7^ 4k, then the solutions (78) appear to be misaligned after i — > i + n. This is not a problem 
because they are on different sheets which are connected by suitable powers of the matrix 
C. When n = 4k we can compare s n+ i with s\. In general, the w variables are related by 
w n+ \ = Wi + Wq. This expresses the fact that there is a logarithmic branch cut. This implies 
that there is a £ dependent monodromy when we relate s^+a ~ s a e~ l a-!i 'o/C-« a «'oC_ if n ^ 4/^ 
then we can choose the origin of the w plane so that there is no ( dependent monodromy. 
Let us now worry about constant parts. These could arise as follows. In the z plane we can 
have a gauge connection whose integral at large z , § A, is non-zero (and proportional to the 
diagonal matrix diag(l, —1, 1, —1)). In writing (77) we have made a gauge transformation 
that sets it to zero. However, this gauge transformation is not globally well defined and it 
appears as an extra gauge transformation that we need to do when we relate s n+ \ with s\. 
When n is odd no such transformation is possible. In fact, it is easy to check that a change 
in normalization of the solutions of the form S{ — > 7 < - _1 ^Sj and S{ — > 7 < - _1 ^ +1 Sj can remove 
any constant from the relation between s n+ i and si. In the case that n is even, there can be 
a constant piece in the formal monodromy. We can take it to have the form s n+ i = 
This is present both for n = 4k and n = 4k + 2. In section 4.5, these properties were derived 
directly from the F-system and the analytic properties. 



D Explicit form of T and ^-functions 

In this appendix we summarize all T and Y functions of the AdS§ problem. We use the shift 
identities (36) to bring the various inner products to the expressions below. In this appendix 
we use the notation (s a Sb) = (s a Sb-iSbSb+i) and (sbS a ) = (sb-iSbSb+iS a ) and drop commas 
inside angle brackets. 
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D.l T-functions 



^2,2fc = ( S -k-2 s -k-lSkSk+l) , T 22 fc+1 = (s_k-2S^k-l s k+l s k+2) 

T\£l = (S-2I-1S2I+1) , T x ^ l+1 = (S-2l-2S2l+l) ) ^l,4«+2 = (■5_2«-2S2«+2) > ^l,4Z+3 = ( ^— 2Z— 3 s 2Z+2 ) 

^3,4Z = (S-2l-lS2l+l) , T 34[+1 = (S-2«-2S2«+l) , 73,4«+2 = (S-2I-2S2I+2) , ^3,4^+3 = ( s -2«-3S2i+2) 

Note the presence of some shifts in the left hand side. Note that Ti jS and T 3i fc only differ by 

Si Si. 



D.2 y-functions 

rr {} _ (S-k-2S-k-lSk+lSk+2) (S-k-lS-kSkSk+l) /Q A\ 

J-2,2k — J-2,2k — JZ r~, Z v , l°4J 

(s_ fe _iS fc+ i) (s_fc_iS fc+ i) 



*2,2fc+l - *2,2fc+l 



(S-k-3S-k-2Sk+lSk+2) {S-k-2S-k-lSkSk+l) 
(S-k-2Sk+l) (S-k-2Sk+l) 



y- _ (S-2l-3S 2 l + 2) {S-2l-2S2l+l} y- _ (g-jj-jggH-l) {S-2l-lS2l) 

1,4 ' +2 (S-21-4S-21-2) (S-21-3S-21-2S21+1S21+2) ' (S-2l-2S-2l-lS2lS2l+l) (S21-1S21+1) 

v (S-2I-3S2I+3) (S-2I-2S2I+2) v \S-2l-lS~2l+l) {S-2I-2S2I+2) 

Y l,4l+3 — 7 Z T 5 r l,4Z+l — 7 T7 T 

\S_2«-4'S-2i-2/ \S-2/-3 s -2/-2'52i+2'S2«+3/ ^2/^2^+2/ \S -21-2S -2l-\S2l+lS2l+2l 

y- _ {S-21-3S21+2) (g— gf— gggf-j-i) ^_ _ (S-2Z-2S2Z+1) (S-2Z-1S2;) 

3,4 ' +2 (s_2«-3S_2«-2S2«+lS2«+2) (s2/S 2i + 2 ) ' 3 ' 4 ' (s _ 2 l- 2 S -2l-lS 2 lS2l+l) (s -21-3S -2l-l) 

v (S-2I-3S2I+3) (S-2l~2S2l+2) v _ \S-2l-lS2l+l) (S-2I-2S2I+2) 

Y 3,4«+3 — 7 T7 = T j Y 3,4«+l 



S-2i-3 s -2i-2S2Z+2 S 2H-3) ( s 2Z+l S 2H-3) (s_2?— 2 S_2Z— 1 ^2^+1 ^2^+2) (S-2Z-3 s -2«-l) 

We also define the functions Y a s to be given by the right hand side of (D.2) without any 
extra shifts. In other words, Y a>s are defined as above with all the products defined at the 
same value of (. In particular, we have included the explicit definitions of Y 2 , s . These are 
related to spacetime cross ratios by replacing — > \ and Sj — > \ iy where Aj are Hodges' 
momentum twistor variables [29], see section 4.7.1. 



E Asymptotic form of the Y functions for the AdS§ case 

In this appendix, we explain how to obtain the asymptotic form of the Y functions for small 
£. For small £ we again diagonalize & z ~ P x ' 4 diag(l, — i, — l,z). We can then use a WKB 
approximation for the study of the solutions. At leading order the solutions are simply given 
by if} a = exp[— i~ a J P(z) 1 ' 4: dz\. We will follow WKB lines where the change of phase is real 
Im(i~ a P(z) 1 ' 4: z/Q = 0. This defines a family of lines. Here it is useful to think about this 
family of lines as lines on the Riemann surface, rather than lines on the z plane. Of course, the 
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sheet > 

sheet 1 > 

sheet 2 ---->---- 

sheet 3 >---- 



(a) 




Figure 11: WKB lines for the AdS§ case. We have four sheets, (a) Shows the conventions 
for naming the WKB lines for each of the four sheets, (b) Through a generic point we have 
four WKB lines passing. The red lines and the blue lines are orthogonal to each other . The 
line for the 1 sheet and the 3 sheet have opposite orientations and are on top of each other, 
here we have separated them for ease in visualization, (c) We evaluate an inner product of 
four solutions by bringing together four of the lines. The order of the lines indicates the 
sign, but we will not keep track of it. (d) When various WKB lines pass through a cut they 
change sheets according to the pattern here, the orientation of the line is preserved. If we 
have several zeros, the effects of this cuts add up. 

Riemann surface is simply x 4 = P(z), and we have the differential xdz = P{z) 1 ^dz. Different 
sheets of the Riemann surface are associated to different solutions of the linear problem. We 
label these sheets as x a , a = 0,1,2,3, in a cyclic way. We also have that x a = i~ a XQ. The 
WKB lines have Im(x a z/() = 0. They are oriented pointing towards the direction in which 
the solution increases. We label the WKB lines for different sheets by different types of 
lines in the z plane, see figure 11(a). Of course, lines for the 1st and 3rd sheet coincide are 
oriented in the opposite way. Similarly for lines and 2. Note, that we also want to keep 
track of the ( independent part which comes from the connection A. The A z component of 
the connection can be diagonalized and set to zero. Only the diagonal components of A z are 
relevant for the small ( WKB approximation. These diagonal components are constrained 
by the Z4 projection condition to be of the form 



The off diagonal components are not relevant for us, since we are neglecting higher order 
terms in the small ( expansion. Again we can think of a z as a one form on the Riemann 
surface. The Hitchin equations imply that da = 0. When we go around a zero of the 
Polynomial P ~ (z — z ) we change sheets figure 11(d). Moreover, going 

around a zero, also a z has a Z 2 branch cut, so that a z — > —a z . In other words, we really 
should think in terms of a U(l) gauge field on the Riemann surface. The Riemann surface 
has a Z4 symmetry x — > ix and a is constrained to be odd under this symmetry. 

The small solutions Si correspond to solutions that are associated to one of the sheets of 
the Riemann surface. In the large z region the sheet of the Riemann surface associated to 
solution Sfc is simply /cmod(4). 

We will draw WKB lines associated to an inner product (sj, Sj, Sk,si) as four lines that 
are incoming if the product is in the numerator of the expression we want to evaluate. Each 
of the lines meeting at the inner produce lives on one of the four sheets, see figure 11(c). The 



A z = a 2 diag(l, -1 



1, — 1) + off diagonal terms 



(85) 
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lines start from each of the asymptotic regions associated to the corresponding Stokes sector 
and they end on the common point where the product is evaluated. We do this with WKB 
lines. Once we identify these lines we can move the point around on the Riemann surface in 
order to simplify the final expression. If the inner product is in the denominator, then we 
reverse the orientation of all the lines, without changing the sheet numbers. 

We take a real polynomial with all its zeros on the real axis and with P(z) > for 
sufficiently large real z. We find it convenient to run all the cuts towards the left of the 
zero of the Polynomial P. The solution s& lives in sheet k (modulo four). It is represented 
as a line coming in from infinity in the corresponding Stokes sector. These lines represent 
integration contours on the Riemann surface where we integrate both P 1 l^dz = xdz as well 
as a. 



>1 



m/4 




Figure 12: Schematic picture of the patterns of WKB lines for a Polynomial whose zeros 
are on the real axis. In (a) we show the pattern when the phase of ( is and in (b) we show 
it when the phase of ( is 7r/4. Lines of different colors are supposed to be crossing at right 
angles. We have five lines ending at each zero. The purple line goes between two different 
zeros. (Its precise type depends on whether we put it above or below the cut). The k th Stokes 
sector corresponds to a WKB line starting on sheet k. Notice that we have only drawn the 
lines that end at zeros. These lines are not useful for evaluating inner products but they are 
lines that separate different families of WKB lines. (The figures in this appendix are best 
viewed on a color monitor). 



With all these preliminaries we are ready to start evaluating some cross ratios. We will 
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evaluate some cross ratios when the phase of £ is 1 and some where the phase of £ = e l7f//4 . 
The patterns of WKB lines for our choice of polynomial for these two cases are given in 
figure 12. We could also alternatively evaluate them all at £ — e t7T ^ 8 where the pattern of 
WKB lines is something intermediate between the patterns in figure 12. 

It is convenient to redefine the Y functions so that they all correspond to some cross ratios 
evaluated at the same value of (. In other words, we define Y asS which are given precisely by 
the expressions in the right hand side of formulas (84) , but with no shifts of the ( argument 
in the left hand side. For instance, Y 2j2 k(C) = *2,2fc(C) and Y 2i2 k+i(0 = Y 22 k+v 



Y 2 ,2 





(b) (c) 



Figure 13: (a) We show the WKB lines that we use to evaluate the cross ratio Y 2 2 . Each 
black dot corresponds to each inner product appearing in the cross ratio Y 2y2 , see (86). The 
inner products that are in the numerator have incoming lines and the ones in the denominator 
have outgoing lines. We have set some dots right on the cuts, we could move them slightly 
up or down, and then we would indeed have four different types of lines ending on them. In 
(b) we have recombined the lines that we had in (a). In (c) we have moved the black dots 
around and we have "annihilated" them. We see the final cycle, 72,2 corresponding to 12,2 • 
Note that the lines going between the two zeros are on the same sheet and have opposite 
orientations. Note also that a would be continuous through the cut going between the 2nd 
and 3rd zero. Thus there is no contribution from a to Y 2j2 . For that reason, there is no 
constant term in (87). 

Let us start with 

Y ~Y - (-3,-2,2,3)(-2,-l,l,2) 
2 ' 2 ~ 2 ' 2 ~ (-3,-2,-l,2)(-2,l,2,3) (86) 

Here we are using a short hand notation (k,l,m,n) = (sk, si, s m , s n ) . The WKB lines 
necessary to evaluate this quantity are displayed in figure 13(a). By reconnecting each pair 
of lines going in and out of an asymptotic region, we arrive at figure 13(b). Then, by bringing 
the four dots together and reconnecting lines on the same sheet and in the same orientation, 
we arrive to the closed contour drawn in figure 13(c). 

The procedure necessary to evaluate any of the Y 2s is basically the same and we get the 
contours displayed in figure 14. To evaluate Y 2>2k+ i we need to evaluate Y 2 2k+X at i(. The 
manipulations are identical, up to a relabeling of the solutions. For ( real we can use the 
WKB lines that are associated to the pattern in 12(b). When we draw the contours in figure 
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(a) 




Figure 14: In (a) and (b) we show the cycles 72,5 that correspond to the various Y 2jS 
functions. We have indicated the sheet where the point with the dot is. Of course there are 
cuts emanating from the zeros, etc. 

14 we also note on which sheet they start. It is enough to label the sheet on which any 
point on the contour is, which determines the rest. It turns out that there is no contribution 
from the gauge connection a for this contour. We can see this by shrinking the little loops 
around the two zeros and then noticing that the two lines that go between the two zeroes in 
figure 14(c) are on the same sheet but have opposite orientation. In addition, we see that a 
is continuous across the cut joining the second and third zero in figure 14(c). This implies 
that the a dependent contribution vanishes. 

In summary, this means that the 1^,^ have the following behavior for small s. 

log y 2 , s = ^ + o(C) Z 2 , s = -£ xdz (87) 

' ^72, s 

where 72 iS are the cycles denoted in figure 14. Notice that the contours all look the same, 
except for the sheet they start on. ( As a first approximation the reader can ignore the 
subtlety about the sheet number where the contour starts.) The sheets change in such 
a way that logY"2,s(C = e 9 ) ~ _ m2 ^ e — with m 2jS real for our choice of the polynomial. 
Here m 2iS are basically the Z 2:S up to possible factors of e 47r//4 (for s odd), see 94. In other 
words, if we consider one of the contours in figure 14, then the four different choices of the 
sheet it starts on would give four possible cycles on the Riemann surface. These would yield 
Z2S, iZ2, s , —%2,s, —iZ2 tS - The asymptotic behavior of Y2, s for ( small and positive is governed 
by the cycle that renders Y 2jS exponentially suppressed. 

We can now evaluate one of the Yi s . Let us evaluate, 

(-2,2)(-l,l) = (-2,1,2,3)(-1,0,1,2) 
1,1 (-2, -1, 1, 2) (0, 1, 2, 3) (-2, -1, 1, 2) (0, 1,2,3) 

The corresponding lines are drawn in figure 15(a). They can be reconnected as in figure 
15(b). And finally lead to the figure eight contour in figure 15(c). It can be seen that there 
is an a dependent contribution to this contour. We can deform the contour to the one in 
figure 15(e). We see that the upper and lower parts of the contours really add, due to the 
properties of a as we cross a cut. Namely, a changes sign across the cut going between the 
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(a) 



(b) 




Figure 15: In (a) we show the WKB lines that are used for evaluating Yi,i- In (b) we 
have reconnected the lines. In (c) we get the final cycle 7^1, which is the figure eight cycle 
indicated. In (d) we see the cycle we would have obtained if we had followed the same 
procedure for Y" 3jl . We see that the cycles for Yi jS and Y 3tS differ by a change in orientation 
and by a shift in sheet by two units. This does not affect the ( dependent piece, but it 
changes the sign of the ( independent piece coming from the gauge connection a. 

first and second zero. This means that the upper and lower contributions of the contour add 
up and give a (generically) non zero answer for a a ^ s . Finally, in figure 16 we draw the cycles 
7i jS for s = 1, . . . , 5. They all have the same eight shape and differ by the sheet, orientation 
and zeros they encircle. 



Yl,5 sheet Yl,3 sheet 2 Yl,l sheet 3 




Figure 16: In (a) and (b) we show the cycles 7i jS that correspond to the various Y 1<s 
functions. We have indicated the sheet where the point with the dot is. Of course, there are 
cuts emanating from the zeros, etc. 

We can similarly evaluate 

Y (-2,2)(-l,l) _ (-3,-2,-l,2)(-2,-l,0,l) 

3,1 (-2, -1, 1, 2) (0, 1, 2, 3) (-2, -1, 1, 2) (0, 1, 2, 3) 
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Figure 17: This figure shows the intersection form for all the cycles associated to the Y 
functions. If an arrow points from Ya to Yb we have (7a, 7b) — 1 otherwise the intersection 
vanishes. 

Notice that this differs from Y\a by an overall bar. This exchanges the top of the diagram in 
figure 16(a) with the bottom. The final contour is shown in figure 15(e). This differs from 
the one for Yi t i in figure 15(c) by the orientation and a change of sheets by two units. This 
combined operation does not change the integral of xdz, which contributes to the l/( term. 
However, it does reverse the sign of the contribution due to the gauge field a. 

Let us make a comment on the result for Y a ^ s as we vary a for a given s. Notice that Yi jS 
is given by a figure eight contour, and so is Y 3tS , but with a different figure eight contour, 
e.g. see figure 15(c) and (d). In principle, the Riemann surface has four figure eight contours 
which start on different sheets. Our course, the sum is zero. In addition, we can have the 
sum of two consecutive figure eight contours. This sum can be deformed into the figure 
"double eight" contour that gives the result for Y 2}S . This implies the pattern of masses that 
we get for Y a ^ s for a given s. 30 

In summary, we get 

log Y a , s = Z f^ + {a-2)c s + o(C) ; c, = - <f a (88) 

We also have that 

Z llia = Z, 3tS , Z l2:s = (1 - (-l)-i)^ (89) 

The last equality follows from the fact that the contours for Y 2jS are a sum of two figure eight 
contours on neighboring sheets. Once we take into account the shifts in ( in the definitions 
of Y a ^ s we see that (89) give the relations between the masses that we found in the F-system. 

It is interesting to compute the intersection form of all the cycles associated to the various 
Y functions. We get the result in figure 17. This result is used in section 3.6. 

30 This is a result of the Z4 symmetry of the Riemann surface. If we were to break that symmetry we 
would no longer have this pattern. 
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The analysis of the large ( behavior is very similar to what we did above. We now 
diagonalize $ s — > P(z) 1//4 (1, i, — 1, i). We have the same pattern of WKB lines and the same 
cycles associated to the Y ajS . But now we will be integrating xdz to obtain Z la s . In addition, 
there is a contribution from the gauge connection A. In this case we can set A z = and we 
focus on the diagonal part of the A s connection, which is given in terms of a one form d z on 
the Remann surface. We have then an expansion of the form 

\ogY^ s = Z^X + (a-2)c s + o(0 

The constant parts c s in (88) and c s are, in principle, different. They are constrained by 
the Y- system equation and we expect that, in the end, we have only one constant per value 
of s which is actually a free parameter. This is most clear for low temperatures (or large 
values of the polynomial P). In this case, we can approximately simultaneously diagonalize 
$ z and $ f . As a result, a and a coincide and c s = c s (see next appendix). As explained in 
the main text, in the general case, the average of these constants is a free parameter while 
the difference is determined by the equations. 

Let us make a comment on the validity of the derivation of the WKB analysis. In this 
case, a given WKB curve gives a good approximation to the cross ratio in a region of size 
7i / 2 in £ around the value of ( for which the WKB curve exists. This is good enough for 
us to derive the F-system, since we need displacements only of tt/4. Note that the WKB 
curves exist for a range of ( of total angle n/2 centered around the values where we did the 
computation. Thus, the total validity of the WKB analysis is a range of ( of angle 7r centered 
on real (, for which we did the analysis. These remarks are valid as long as the polynomial 
has its zeros along the real axis. 

When we move the zeros away from the real axis, the quantities Z a s become more general 
complex numbers, though still constrained as in (89). Thus we have 2(n — 5) complex 
parameters which correspond to the motion of the zeros of the polynomial. (There are other 
(n — 5) parameters in the constants c s + c s ). Now, as we move the zeros of the polynomial 
we will find that at some point some of the WKB lines will cease to exist. However, having 
derived the integral equation in one region in parameter space (i.e. for real masses), we can 
move to other regions by analytic continuation. In doing so, we have to change the integral 
equation according to an easily derived pattern, as explained in appendix B. 

F Reality conditions for the Y functions 

The reality conditions for the Y functions depend on the signature in which the solutions 
are embedded. More precisely, for SO(p, 6 — p), they depend on the parity of p. For p even, 
which corresponds to a boundary of AdS with (3, 1) or (1, 3) signature, the reality conditions 
in a suitable choice of gauge have the form 

4 = -BA Z B' 1 , $t = BQzB- 1 , (A(C)) j = -BAi-hB- 1 (90) 
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For a suitable B for each case. For p odd, which for instance corresponds to the case of (2, 2) 
signature, the reality conditions have instead the form 

At = BA z B~ l , = B$ z B~ l , (.4(0)* = BA^B' 1 (91) 

For a suitable B for each case. The difference between different parities of p comes from the 
fact that conjugation interchanges the spinorial and anti-spinorial representations for p even, 
but it doesn't for p odd. In addition, we also have the Z 4 projection condition 

A(i() = -CA^QC- 1 , A(-C) = U- l A{C)U , U = C t C- 1 (92) 

where the second equality arises by applying the first one twice. 

Let us first focus in the case in which p is even. If V'(C) * s a solution to (d + A(())ip = 0, 
then by complex conjugating the above equation and using (90) (92) we get that 

(d - ^(i/c*) w(or = o , d = c-wb* 

This implies that D(if)(£))* oc If we take if) to be a small solution Sk, then we see 

that (sfc(C))* oc D _1 Sfc(l/C*)- Note that D will drop out when we compute inner products, 
and any constant of proportionality drops out when we compute cross ratios. Thus the effect 
of performing a conjugation is to replace Sk by Sk in all formulas. This implies that 

(Y ajS (or = >wi/a 

In particular, for large 9 = log ( we see that 

[log(y ljS /y 3)S ) (+oo)]* = -log(Y 1>8 /Y 3t9 ) (-oo) 

which implies that c s = — (c s )*, or C s purely imaginary and D s purely real. 

For p odd the situation is different. By complex conjugating (d + A(())ip(() = 0, and 
using (91), we simply obtain (ip(())* oc if) (!/(?). This implies (sjfc(C))* oc s^(l/C*) and hence 

(Y a>s (or = Y ajS (i/o 

implying the Y's are real when ( is a phase. In this case c* — c s , and C s is real and D s is 
purely imaginary. 



G Components of the full area 

As seen in the body of the paper, in order to compute scattering amplitudes at strong 
couplings, we need to compute the area of minimal surfaces in AdS, given by 



A = J d 2 zTr 
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Since for solutions relevant to scattering amplitudes, Tr [Q z <& z ] ~ (PP) 1 / 4 , this area diverges 
and needs to be regularized. This can be conveniently done by dividing the area in different 
contributions (we refer the reader to [17, 15] for the details). Below we give the results for 
the case in which n ^ 4k, where the treatment is simpler 31 

A reg = J d 2 z (Tr - 4(PP) 1 / 4 ) 

A periods = 4 f d 2 z(PP) 1/4 -A j d 2 w = 4 f d 2 w-4 f d 

J J So J S J So 



2 w 



^cuto ff = 4 / d 2 W 
' °~0,ZAdS>£ 



where we have defined the w— plane by dw = P(z) l l 4 dz. A reg is the non-trivial function 
that is computed by the free energy of the TBA equations and is the central quantity of this 
paper. 

In A periods , £ denotes the surface defined by the polynomial P(z), while S is a reference 
surface with a single branch point at the origin and the same structure at infinity. When 
computing it, it is important to use a cut-off in the if— plane, such as \w\ < A, with A ^> 1. 
It can be evaluated in terms of the periods of the Riemann surface x 4 = P(z) and we give 
its explicit expression below. 

When computing A cut - ff we impose a physical cut-off, namely that the radial AdS 
coordinate should be larger that certain small e. It can be conveniently written as the sum 
of two contributions 

■A-cut-off — g ^ l°g 2 ( e2x 4 2 ,i+2) + A B DS-like 

i 

the first term is the standard divergent term, expected for light-like Wilson loops/amplitudes 
[43, 9]. The second term is 



.4 



BDS-like 



1 n / 2K \ 

- ^ log 2 xJ J+2 + ^(-l) fc+1 logx 2 i+2 logx 2 +2fe+M+2fc+3 , n = AK + 2 



k=0 
2K 



= log 2 x 2 jl+2 + f](-l) fc+1 logx 2 ii+2 logx 2 +2fe+M+2fc+ 3 , n = AK + 2±l 

t=\ \ fc=0 / 

AsDS-uke is a finite term which obeys the conformal Ward identities of broken conformal 
invariance [10] . Actually, it is the unique solution of the anomalous Ward identities that 
can be written only in terms of next to nearest distances x 2 i+2 . 

It is customary to subtract the one loop result Abds written down in [9], with the 
appropriate overall factor. As both, AsDS-iike and Abds satisfy the Ward identities, their 



31 The case n — 4fc is subtle, since the w— plane possesses a monodromy at infinity. This was explicitly 
treated in [17] for AdS^ kinematics. We expect similar results for the AdS§ case, but we have not worked 
out the details. It would be interesting to work them out for n = 4fc in AdS$. 
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difference is a function of the cross-ratios. In practise, in order to express A BDS — A BDS -uke in 
terms of cross-ratios, one simply starts with any non-next to nearest distance x? ■ appearing 
in Abds an d write it in terms of the unique cross-ratio which involves x|- and next to 
nearest distances x^ i+2 . The strong coupling answer, however, organizes more naturally as 
described above. 



G.l Expression for A perio( i s 

As already mentioned A perioi i s can be evaluated in terms of the periods of the Riemann 
surface x 4 = P(z). More precisely, it is given by 



A 



periods 



(93) 



where 7 denotes the collective pair (a, s) and tf 7i7 ' is the inverse of the intersection form of 
the cycles. Such inverse exists for the case in which n is odd. Using the relation between the 
different periods and the relation between these and the masses 

Z hs = Z 3 , s = Z„ Z 2 , s = (1 - {-l) s i)Z s 

m 2s +i = -2Z 2s+1 , m 2s = -2e~^Z 2s (94) 

we can write A per i d s in terms of the masses of the TBA equations (47). We obtain 

Aperiods K^ijTfliTflj 

The form of /C depends on the parity of (n — l)/2. For instance, for the case n = Ak + 5, K. 
is a matrix of 4k by 4fc and is equal to 



1 

71 



(0 

1 


1 
2 
1 


1 
1 


"A 
1 


2 


v/2 
1 
1 

v/2 


1 


f 


V 2 


f 
2 


2 



:J - (Xi)w = e(j-i)(-i) < - J ' +1 l e(o) = o 

/O 1 >/2 \ 

^3 1 



(95) 



2v/2 



0^2^ 
1 V2 1 
\0 / 



The case n = 4/c + 3 has very similar expressions if we write /C in terms of a 4fc by 4k matrix 
whose last two rows and columns have to be chopped off. The expression coincides with that 
in (95), with different matrices K.± and /C3 



ICi 



2V2 



fy/2 1 -1 \ 

1 V2 -y/2 

-1 

\-l -y/2 -1 / 



/c. 



2V2 



fy/2 1 -1 \ 

2 V2 -V2 

>/2 1 -1 

\ 00 / 



The case n = 4k + 2 can be regarded as a limit of the case n = 4k + 3, in which we take one of 
the zeroes of the polynomial very far away. The correct prescription is to start with Apl^ ds , 
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subtract ^(m 4fc „ 2 — ^{jn^k-2 — and then take |m4fc_ 2 | very large, k is then uniquely fixed 
by the requirement that linear divergences cancel. Again, the result can be written in exactly 
the form (95), with different matrices and then chopping away the last three rows and 
columns 



/Ci 



1 

4 



/-l 1 V2\ 

V2 2 

1 V2 1 V2 
\V2 2 V2 / 



/C, 





/ 1 





-1 


-V2\ 


1 







-V2 


-2 


4 


1 





-1 






V o 








o / 



In particular, it can be checked that we reproduce the correct answer for the hexagon case 
[15]. 

The result for A peTio( i s is particularly simple in the case of AdS%. It can be computed by 
using Ap^ ds = —iw rs Z r 2 s , but in this case the inverse matrix w r ^ s is much simpler. For 
n = 2h, gluons, with h odd, and defining m<ik = rri2k, fh\ = rti\ and fhik+x = ^2fc-i + 
for k — 1, we obtain 



h — 3 
2 



1 

k=l 



(96) 



where we have used the relation between the periods Z r and the masses m r from section 3. 



H Direct computation of the regularized area 

In this appendix we consider the regularized area A reg for a particular class of regular poly- 
gons which can be embedded either in AdS 3 or AdS^ and correspond to special radially 
symmetric solutions (in a sense which will be clear momentarily). We will then compare 
such results with the answer for the free energy of the Y-system in the high temperature 
limit. 

We will be interested in the case of (2, 2) signature, since in this case we can embed 
both kind of solutions. Strings on AdS^ can be described in terms of the usual holomorphic 
function P(z) and two fields a and 0. We can choose a gauge in which the connection 
becomes 

A _ 1 (-da - d?a 3 \ _ _ _ / e-^V2P(z)^ L%, <~™)\ 

z 4 V da- df3a 3 J ' ^ z " \^ J 

And A 2 = —Al, $ s = $t We have written the connection in terms of two by two blocks. We 
will consider a symmetric configuration in which all the zeroes of P{z) are together, namely 
the holomorphic function is a homogeneous polynomial P(z) ~ z n ~ 4 . In this limit, all its 
periods, and hence the masses entering the TBA equations, vanish. On the other hand, for 
this particular case, it is consistent with the equations of motion and boundary conditions 
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(Xi ,x 2 ) 



5 6 

O 



Figure 18: Two classes of regular polygons which live in (2, 2) signature and can be embedded 
in AdS 3 (a) and AdS^ (b). Both posses a Z n symmetry. 

to set a and (3 to be functions of the radial coordinate only, \z\ or \w\ = p. It is convenient 
to write the equations they satisfy in the w— plane, defined by dw = P(z) 1 ^ 4 dz 

a"(p) + _ 8e & + 8e~ & cosh (3 = 

P 

/3"(p) + — -8e-"sinh/3 = 
P 

where we have defined the shifted field a = a — | log |-P(^)| — log 2. For the scattering of 
n gluons we will consider two different configurations. The first configuration, which exists 
only for n even, corresponds to regular polygons that can be embedded into AdSs and hence 
(3 = for these. Geometrically, the span a regular polygon of n sides in the (xi,x 2 ) plane 
and a segment in the (ti,t 2 ) plane, see figure 18 (a). They have been considered [17] and 
their area has already been computed there. 

The second configuration corresponds to regular polygons that can be embedded into 
AdS±. Geometrically, they have the shape of a regular polygon of n sides in the (xi,^) 
plane and a regular polygon in the (ti,^) plane, of n/2 sides if n is even, or n sides if n is 
odd, see figure 18 (b). As already mentioned, we expect a and (3 to be radially symmetric. 
Furthermore, following the discussion in appendix B of [15], we expect the following boundary 
conditions at the origin 

n — 4 

a = 2 logp + c a + ... 

n 

n — 4 

(3 = 4 logp + c/j + ... 

n 

and a exponential decay at infinity. Finally, we are interested in the regularized area: 

A reg = Inn \ pdp(e a — 1) 



(ti,t 2 ) 

1,3,5 

(a) 

2,4,6 
1,4 

^>2,5 (b) 
3,6 
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It is at present unknown how to compute analytically the above area. Based on our experience 
with the first class of regular polygons, we expect a factor of tt, times a simple rational 
function of n. Furthermore, this expression should vanish for n = 4 and, in addition, we also 
know the correct value for n = 6 from the results of [15]. 

We have solved the above equations numerically. The constants c a and cp are unknown 
and are fixed by the constraint that the solutions decay at infinity. One should solve a 
shooting problem, trying several values for these constants approaching the correct ones by 
requiring that the solutions decay exponentially at infinity. We performed such procedure 
for n = 7, 8 and n = 10 and obtained the following numerical results 

A™=J « 4.882, A?= 8 « 7.067, A?= 10 « 11.85 

Note that these expressions are well approximated by ||7r, |7T and ^-tc. Once the regular- 
ized area is computed, in order to compare it to the free energy, we need to subtract an 
n— dependent constant, ensuring that the free energy vanishes for the situation in which all 
the zeroes are well separated. In this case, the regularized area approaches the quantity of 
zeroes, i.e. n — 4, times the contribution of the regular pentagon, which corresponds to a 
single zero. 

Extracting from [15] the value of the regularized area for the regular pentagon, A penta = 
|7r, we can give the final formula for the direct computation of the free energy 

3 

A free = A reg - (n-4)-7r (97) 

o 

From the numerical computations we can guess a very simple final formula for the regular 
polygons of the second class: 

l(2 ) l(n-4)(n-5) , . 

41=2 ^ (98) 

This vanishes for n = 4, as expected, gives the correct result for n = 6 and agrees very well 
with the numerics for n = 7, 8, 10. 

Let us now plug in (97) the results for the regular polygons that can be embedded into 
AdSs, read off from [17]. Substracting the appropriate contribution we obtain the appropriate 
quantity to be compared with the free energy for regular polygons of the first class: 

~^~free ^-sinh ij^ 4) Ap en t a g 0n 7T — (99) 

where A Sinh is the result in [17]. We see that (99) vanishes for n = 4, as expected, and also 
agrees with the answer of [15] for n = 6. 

AdS 3 limit 

In the AdSz limit the expression for the regularized area of regular polygons is of course 
the same as for the polygons we were discussing above (which could be embedded AdSs). 
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However, in order to make a direct comparison with the free energy we need to substract a 
different contribution. Considering n = 2n gluons, the holomorphic polynomial has n — 2 
zeroes, each of which gives a contribution equal to that of the regular hexagon, we obtain: 

4K? = >W-(ft-2)^r = ir <"-W B - 4 > (100) 

Note that (100) is different from what we would obtain from the free energy in the AdS§ 
case (99), due to the fact that we subtract different contributions. 



I Regular polygons 

In the body of the paper we have solved the TBA equations in the high-temperature limit, 
in which constant solutions of the Y-system are relevant. In this limit, and for an even 
number of sides, there is a family of solutions parametrized by a single parameter \x. The 
solution interpolates between the two kind of regular polygons described in appendix H. In 
the following we describe the geometrical picture of such polygons. As we will see, this allows 
to compute analytic expressions for the various cross-ratios. These expressions can then be 
compared to the results obtained from the Y-system equations testing both, the geometrical 
picture and the Y-system equations. 

In addition, when the number of sides is odd, there is a discrete family of regular polygons 
that can be embedded into the boundary of AdS±. We present them in the second subsection 
of this appendix, and compare their cross-ratios with the respective solutions from the Y- 
system. Finding again agreement. 



Regular polygons with an even number of sides 

In appendix H we have seen that there are two families of regular polygons with an even 
number of sides, which can be embedded into the boundary of AdS% and the boundary of 
AdSi respectively. It is possible to construct a family of polygons that interpolates between 
these two in the boundary of AdS$. For that it is convenient to describe the boundary of 
AdS$ in terms of projective coordinates with (2,4) signature 32 . 

-Z X Z X - Z 2 Z 2 + Z Z Z Z = (101) 

We propose that the location of the cusps of the polygons under consideration is 

Z(=(-iy£ 1+ i, Z» = £ 2 e^, Zl = i z e^ (102) 

where £1,2,3 have to be chosen in such a way that (101) is obeyed and the distance between 
two consecutive points is light-like, namely 

1 + £ 2 - £ 2 + £ 2 = , -t\ + sin 2 f— } t\ - sin 2 (-) t\ = (103) 



:i:2 



These polygons live in a boundary with (1, 3) signature. 
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This leave us with one free parameter which we can parametrize in terms of an angle as 




where <f> runs from zero to ^i^7r. This parameterization is engineered in order to ease the 
comparison with the results in the main text. As explained below the formal monodromy \x 
appearing in the main text is given by \x = . 

We can interpret the solution as in (1,3) signature, with the points living on a spatial 
S 1 times a temporal S 3 . When <p = 0, £ 2 takes its maximal value and £% = 0. All the points 
on the S 2 C S 3 are located on its equator, forming a regular polygon of n/2 sides, and we 
recover the regular polygons embedded into the boundary of AdS^. When <f> = ^§-71", £ 2 
vanishes and £\ takes its maximal value, so the points in the S 2 C S 3 alternate between 
the south and north-pole and the solution reduces to the usual regular polygon that can be 
embedded into the boundary of AdS 3 . For intermediate values of 0, we interpolate between 
the two solutions and the corresponding minimal surface span the full AdS$ 33 . Furthermore, 
note that the solution possesses a Z n symmetry. 34 

In order to compute the cross-ratios corresponding to this solution, and the relation 
between \i used in the body of the paper and the parametrization used here, it is conve- 
nient to compute the spinors \ corresponding to this solution. In other words, we need to 
give a description of the solution in terms of twistors. The space-time coordinates can be 
parametrised in the spinorial representation as follows 

/ z x z_ 2 z 3 \ 

-Zx Z 3 Z 2 
-Z 2 -Zj -Zi 

\-z 3 -z 2 z x / 

we will denote by a, = 1, ...,4 the spinorial indices. The points X 1 satisfy the following 
conditions 

e a ^ s X^X; s = , e a ^ s Xl X\f = 
As a consequence, they can be written in terms of spinors A^,, such that 

~yi _ \i— 1 \i \i \i— 1 
— A a A /3 ~ A a A j3 

See figure (19). There is a simple recipe to determine the spinor A 1 , since we know it has to 
satisfy the following equations 

e «/*y*X* 7 AS = , e a ^ s Xi+ l \i = 

given the points X % and X l+l , each of the equations above specifies a two dimensional plane 
and the intersection of these planes give us A*. Hence A* will be completely specified by 

33 The full ^ciS^ is spanned since the radius of the temporal S 2 changes as we go between the two endpoints 
of the interpolation. 

34 More precisely, two consecutive points are related by an 0(6) (but not 5*0(6)) rotation X t+1 = RX l , 
where R does not depend on i. 
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Figure 19: We can associate a twistor A* to each edge of a null polygon, such that the position 
of the cusps is given by \ l ~ 1 \ t 

the above equations up to an overall normalization factor. Of course, these factors are 
inessential, as they will drop out from the computation of any cross-ratio. We obtain the 
relatively simple result 

A oc {«i (1 - a k ) (1 - b k )e « , n 2 e « , k 3 {1 - a k )e » , (1 - b k )e »} 

where 

■ 27T k 7T (j) 

a k = i — tan — tan—, b k = (— 1) tan — tan — , 
n n n n 

and Ki are some fc-independent constants. The spinors posses a Z n / 2 symmetry \ l+2 oc UX l 

with [/ = diag(e~ : ). Having the spinors we can compute the invariants 

(i,jXl) = ^ 5 X a \^\% (104) 

In order to compute the invariants entering in the cross-ratios Y s ^ m , it is sometimes convenient 
to have A^ = e^Ajp^A^ 1 . We obtain 

— 7, 3t7rk . \ / \ 3a7rfc . j \ ink . . iirk -~ 

A oc {/tie « , K2 (1 + a fe ) (1 + b k )e « , k 3 (1 + fc )e •> , (1 + a fc )e « } 

For some constants fiij. Finally, using the explicit expressions for the cross-ratios Y S)Tn in 
terms of invariants, we can obtain analytic expressions for them. These are ratios of four 
such terms such that any Aj (or Aj) will appear as many times in the denominator as in the 
numerator. In particular, such cross-ratios do not depend on the k independent constants 
Ki,Ri. In other words, these constants can be set to one by a conformal transformation. 
Re-shuffling the elements of the resulting spinors Aj, we see that they precisely agree with 
the ones used to write the solution in the high temperature limit of the AdS^ Y-system, 
equation (64). 

In general the expressions for the V-functions are not particularly illuminating. In some 
cases there is significant simplification. For instance, for the case n = 12, it greatly simplifies 
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and we obtain 

^ 3 = 12 = ^ + ^(3 + V3) e -^cos^ 

Before proceeding, let us note the following fact. Taking — » — <p effectively interchanges 
A -H- A. From the geometrical point of view, this simply amounts to take Z\ — > —Z\ (keeping 
2^2,3 fixed ) and is a symmetry of our configuration. From the point of view of the Y-system, 
it interchanges Yi jS -B- Y 3 s and takes \i — > fi^ 1 , which is also a symmetry of the equations. 

Regular polygons with an odd number of sides 

For an odd number of sides there exist discrete families of regular polygons that can be 
embedded into the boundary of AdS<±. These polygons can be conveniently written using 
(2, 2) signature and are parametrized by the number of sides n and an extra parameter 
r = 2, (n — l)/2. The cusps are located at 

_ 2-Kir 2-ki 

Zi = % , Z 2 = £ 2 e n , Z-i = £ 3 e « 

£2 and £3 have to be fixed in such a way that the distance between consecutive cusps is 
light-like, and to obey (101). 

One can proceed along the lines of the previous subsection, going to projective coordinates 
and computing the spinors. The expression for the spinors are not particularly illuminating, 
however, they satisfy the following relation 

A + = diag(e ™ ,e « ,e ™ , 1)A 

up to arbitrary constants. This relation, allows us to write all the spinors in terms A , which 
can be set to A = (1, 1, 1, 1) by conformal transformations (this is analogous to setting 
Ki — 1 in the previous subsection). In addition, we multiply each spinor by an overall phase 
e -i7r(r-i)-_ ^y" e arr j ve t the very simple result 

)k = ^ e in(r+l)^ e -in(r+l)^ e i7r(r-l)^ e -i7r(r-l)^ 

Furthermore, one finds that A and A agree (up to a multiplication by a constant matrix). 
The computation of cross-ratios is now straightforward and one can check that these are 
indeed solutions of the Y-system equations. 
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